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Preface 

As I always believe that the course about equations from mathemat- 
ical physics is of fundamental importance to every student in college. 
However, it will be a hard time for us to firstly learn it. Can you image 
what is up when you enter a new world where you know just a little? 

I want to point out some key points in our course at this very be- 
ginning. To us, solving a problem means that we reduce the problem 
into some easy problems which can be solved by us. For solving our 
problems with partial differential equations (in short, PDE), we reduce 
them into ordinary differential equations (in short, ODE), which can 
be solved by us. 

We shall derive some typical partial differential equations from the 
conservation laws of mass or energy or momentum. Such equations are 
called equations from mathematical physics. The basic principle used 
here is that the change (rate) of the physical quantity in a body in 
space is created from the change on the boundary and by the source in 
the interior of the body. Hence, we need to use the divergence theorem 
from the differential calculus. 

The equations from mathematical physics are very important objects 
in both mathematics and applied mathematics. However, many of them 
are nonlinear equations. Even for linear partial differential equations, 
we can not expect to have a unifying theory to handle them. Therefore, 
we have to classify the equations and treat them case by case. The 
crucial observation is that many equations enjoy symmetry properties 
like translation invariant or rotation invariant, which will lead to many 
identities for solutions. 

The object in our course is to try to explain some elementary part 
of solving simple linear partial differential equations of first order and 
second order. The principle for studying these equations is reduce them 
into linear ordinary differential equations, which can be solved by tools 
from the differential calculus and linear algebra. 

The most beautiful method for solving linear partial differential equa- 
tions or systems of first order is the characteristic line method. 

The most powerful methods for solving linear partial differential 
equations or systems of second order are the methods of separation 
for variables (also called Fourier method ) and Fourier transformation 
method, which can be used to translated the problems into ODE again. 
Theory of Fourier transformations is one of most beautiful mathemat- 
ical theories so that there are many books about Fourier transforma- 
tions. The other method is to construct the fundamental function to 
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the equation with the homogeneous boundary condition and initial con- 
dition. In the latter method, the divergence theorem plays an essential 
role. 

There is another important tool in the study of heat equation and 
elliptic equations. It is the maximum principle trick with its brother 
called the comparison lemma. This trick is also very useful in the study 
of nonlinear parabolic and elliptic equations or systems. People usually 
use the maximum principle to derive gradient estimates or Harnack 
inequalities of non-negative solutions, which lead to Bernstein type 
theorem. 

To understand the behavior of the solutions to PDE, we use diver- 
gence theorem, Fourier method, or the fundamental solutions to set up 
energy bounds for solutions. The nice forms of energy bounds for so- 
lutions are monotonicity formulae or Pohozaev type identities , which 
can be used to get Liouville type theorem or regularity theorem. 

Energy bounds or gradient estimates are important to understand 
the compactness of solution spaces. 

The aim of this course is to let our students/readers know the basic 
ingredients in the study of simple linear PDE's of first order and second 
order. Sections with * ahead should be skipped by first reading. Just 
as in the study of other courses in mathematics, one can only get 
a good understanding of our subject by doing more exercises 
and by reading more related books. 

I hope you understand my words here. What did I feel in my college 
life? I did feel so lucky that I had attended many courses in math- 
ematics when I was in Nankai University. This experience makes me 
can understand many deep lectures from some famous guys. 

In conclusion, for our beginners, I hope you remember that we try to 
understand PDE by ODE! In other words, before we do "PDE", let 
us do its partner "ODE". Do you remember an old Chinese proverb? 
Here it is, "/ suddenly find that I know the world after I hear your 
words" . 



Written by 
Li Ma 
in 5-1-2007 



EQUATIONS FROM MATHEMATICAL PHYSICS 



7 



Notations and conventions: 

The functions used in our course are bounded smooth function unless 
we point out explicitly. If the domain is unbounded, we assume that the 
functions is decay to zero at infinity. In other word, we assume the nice 
conditions for functions such that we can perform the differentiation 
into integrands or the change the order of integration of variables. 

We denote by R™ the Euclidean space of dimension n. Let x = 
(xi, ...x n ) the standard coordinates in R n . 

We denote by C 2 (D) the space of twice differentiable smooth func- 
tions inside the domain D~R n . 

We denote by C£°(D) the space of infinitely differentiable smooth 
functions with support inside the domain DH n . Here the support of a 
function / is the closure of the set {x G D; f(x) ^ 0}. 

V/: {<),,/ i),f). 

v,/- d xj f. 

3 

is the radial derivative of the function /, where r = \x — p\ is the radial 
coordinate with respect to the center p = (p>) e R™. 
We denote by 

di k) f = k = h + ... + k n 

the k-th derivatives of the function / with respect to the variables 

X (Xi , . . .Xn) . 

PDE=Partial Differential Equation. 
ODE=Ordinary Differential Equation. 
F-transform=Fourier transform. 
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1. LECTURE ONE 

1.1. Divergence theorem. 

Let (x % ) be the standard coordinates in the Euclidean space W 1 . We 
recall the classical divergence theorem: 

Theorem 1. Let D be a bounded domain with finite piece-wise C 1 
boundary 3D in W 1 . Let X : D — > W 1 be a smooth vector field on the 
closure D of the domain D. Let v be the unit outer normal to dD. 
Then we have the following formula: 



I divXdx — X ■ uda. 

J D JdD 



Here divX = £ . f£ for X = (X*) , v = (u t ) and X ■ v = £.XV 4 . 
If X % = J^- for some smooth function, i.e., X is the gradient Vw : = 
(d x iu, ...,d x nu) of the function u, then 

divX = = Au 

*-r* [dx 1 ) 2 

which is the Laplacian of the function u. 

A special case we often use is that X has only one non-zero compo- 
nent, say j — th component, 

X = (0,...0,m,0,...0) 

for some smooth function u G C 1 (-D) fl C(D). Then, divX = dju and 



/ djudx = / uvjda. 

JD JdD 



The classical divergence theorem is a basic theorem used in our 
course. 

We ask the readers to remember this basic formula. Note that when 
n — 1, the divergence theorem is the famous Newtonian-Leibniz 
formula: 

f(b) - f(a) = [ f\x)dx. 

J a 

In dimension two, the divergence theorem is Green's theorem in dif- 
ferential calculus. Let D be a bounded smooth domain in the plane 
R 2 . Recall that the Green formula is the following: for any smooth 
function P, Q on the closure of the domain D, we have 



/ (Q x -P y )dx= / Pdx + Qdy. 

JD JdD 
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It is clear that the form of divergence theorem in dimension two 
is different from that of the Green formula. We want to prove that 
they are equivalent each other. Assume that the boundary dD is given 
by the smooth curve p{s) = (x(s),y(s)), where s is the counter-clock 
arc-length parameter. Write by x' = dx/ds, etc. Then T = (x',y') is 
the unit tangent vector along p(s). Then N = (—y',x') is the interior 
unit normal to the curve p(s) such that {T, N} forms a right-hand 
orthonormal basis system. The outer unit normal to the boundary dD 
is v — —N = (y', —x'). Choose X 1 = Q and X 2 = —P. Then we have 



divX — Q x — P% 



and along the boundary dD, 

Pdx + Qdy = (Px' + Qy')ds = X ■ vds 
Hence, Green's formula implies that 




divXdx = 




which is the divergence theorem in dimension two. 
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1.2. One dimensional wave equation. 

We consider a vibrating string of length I with its two ends fixed at 
x = and x — I. Let u = u(x, t) be the position function of the string 
at (x, t) in the direction of y-axis. Assume that p = constant is the 
mass density of the string. Then the kinetic (or temporal) energy of 
the string is 



1 



T ( u ) = 2 / P u t dx - 



The potential (or spatial) energy is 



pu x dx, 

where c is a physical constant of the string. 

Assume that the string is acted by the exterior force with force den- 
sity fo(x,t) at (x,t). 

Then the energy obtained by the exterior force is 

mdx. 



U(f) = - I hi 
Jo 

Then the Hamiltonian action of the string is 

H{u) = [ dr(T(u)-U(u)-U(f )) 
Jo 



'0 

- [ dr [ pu+dx 
2 Jo Jo * 



t 

2 







dr / pu x dx + dr foudx. 



The Hamilton principle implies that the motion of the vibrating string 
obeys 

8H(u) = 0, 

that is, for every f = £(x,t) E C 2 ((0,/) x (0,*i)), 

<5/J( M ),0>:= ^-H(u + eOU=o = 0. 
de 

We compute that 

<5H(u),(f)> = dr (pu t (f) t 
Jo Jo 
- cpu x <p x + fo<p)dx 



dr (-pu tt + cpu xx + fo)4>dx 
Jo Jo 
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Using the variational lemma (see Appendix), we then get 

u a - cu xx = fo/p, 

which is called the one dimensional wave equation. 

By physical testings, we can find the initial position and velocity of 
the string. So we know that 

u(x, 0) = <f>(x), 

and 

u t (x,0) = ip(x); 
and the boundary condition 

u (0,t) = = u(l,t). 

The higher dimensional wave equation in its simple form is 

u tt (x, t) = a 2 Au(x, t) + f(x, t), 

which is a typical hyperbolic partial differential equation of second 
order. 

We can also add other type of initial data and boundary data for the 
wave equations. 
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1.3. Conservation of mass. 

Let p = p(x, t) be the mass density of a moving fluid in the region 
W. Let u be the velocity field of the fluid. Let v be the outer unit 
normal to the boundary of the domain W. 

Take a small part dx in W. Then the mass in dx is pdx. Hence, the 
total mass in W at time t is the integration of pdx over W, which is 



pdx. 

>w 

Then, its change rate at time t is 

d 
dt 



pdx. 

w 



Assume that there is no fluid generating source. Then the change 
rate is only across fluid rate from the boundary. Take a small part da 
in the boundary dW. Then the amount of — pu ■ uda flow into W. So 
the total amount flowing into W across the boundary is the integration 
over dW: 

— pu - vda. 

JdW 

So we have 

/ pdx — — pu ■ vda. 

'W JdW 



dt 



Using the divergence theorem, we have 
f d 

/ (— p + div(pu))dx — 0. 
Jw 

Since we can take W be any sub-domain, we have 

d 

—p + div(pu) = 0. 

This is the continuity equation in Fluid mechanics. 
If u = (a(x 1 ), 0, 0), then the continuity equation is 

Pt + d x (a(x)p(x,t)) = 0, 

where x = x 1 . 
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1.4. -k What are important for our equations? 

What I want to show here is that once we have a pde from scientific 
background, we must try to solve it and get the explicit expression for 
the solutions in terms of initial data and boundary data. We also need 
to know the uniqueness and stability of our solutions to our equations 
or problems. 

Here we may introduce the concepts of Hilbert space and normed 
space (see any textbook about Functional Analysis). A norm | • | on a 
vector space E is a function from E to non-negative numbers R + such 
that 

(1) 

1/1=0 

if and only if / = 0; 

(2) \af\ = \a\\f\ for any real number a; 
(3) 

!/ + <?!< I/I + l<?l- 

A vector space E equipped with a norm | • | is called a normed space. 
If a sequence (fj) C E satisfy \ fj — fk\ — > as j, k — > oo, then we call 
(fj) a Cauchy sequence in (E, \ ■ |). If every Cauchy sequence has a 
limit in E, then we call (E, \ ■ |) a Banach space. If further, the norm 
| • | comes from an inner product < •, •>, that is, |/| 2 =< /, / >, we 
call (E, <•,•>) a Hilbert space; for example, 

<f,9>= / f(x)g{x)dx, f,geL 2 (D), 

J D 

By definition, we say that u(x, t) is the stable solution (in the norm 
| • |) to the problem 

A[u](x,t) := (d t u- Lu)(x,t) = 0, 

where t G [0,T], and 

k 

L:=P(x,d x ,...,di k >) = J2*Ax)d ( J ) 

3=0 

is a partial differential operator of the variables (x), with initial data 
u(x, 0) if for any e > 0, there exists 5 > such that for all solutions 
w(x,t) to A[w] = with initial data w(0), 

\w(0) - u(0)| < 5, 

we have 

sup \ w(t) — u(t)\ < e. 

t 



14 



LI MA 



There are also many other definitions for stabilities, which depend on 
what we wanted in reality. 

To get uniqueness and stability of solutions to evolution problems, 
we need to find some nice norm for solutions and derive a Gronwall 
type inequality. For stationary problems, we have to use maximum 
principle trick or monotonicity formulae for solutions. 

In our course, we need to know some typical methods to solve linear 
pde's of first order and second order. Historically, the solution expres- 
sions were guessed by scientific workers from physical intuition. So it 
is important for us know the scientific background of the equations. 
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2. LECTURE TWO 

2.1. Heat equation and Laplace equation. 

Heat equation and Laplace equation are very important objects in 
PDE theory. 

We use the energy conservation law to derive the heat equations. 
Let u be the heat density of the body Q, and let q be the heat flowing 
velocity. By Fourier's law, we have 

q = —kVu, 

where k > is a coefficient of physical quantity. 

Let fo(x,t) be the density of heat source. Then as in the case of 
continuity equation, we have, for a physical constant c > of the 
region Dell, 

d t / cpudx — q ■ (-u)da + / / pfodx. 

JD JdD J JD 

Using the divergence theorem we have 

/ Q ' {—v)d<J = — div(q)dx. 

JdD JD 

Hence, since D can be any domain in the body, we have 

cpd t u = —divq + pf . 

Using the Fourier law, we have 

cpd t u = divikVu) + pf , 

which is called the heat equation with source f . Assume that c = 1 = 
k = p and / = 0. Then we have the heat equation 

u t = Au, (x, t) GO x (0, oo), 

which is a typical parabolic partial differential equation of second order. 

Physically, we can measure the heat density on the boundary dQ, 
and we may assume that u(x,t) = for (x,t) G dQ x (0, oo). The 
initial heat density at t — can also be detected. So we assume that 

u(x, 0) = <f>(x), x e Q. 

The heat equation coupling with initial and boundary conditions is 
called a determined problem of heat equation. 

Of course, we may have other kind of boundary conditions like 

du 

k—(x,t) — au(x,t), (x, t) G dVL x (0, oo). 
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When we have stationary heat in f2, that u t = 0, we have the follow- 
ing equation: 

Au + / = 0, 
where / = pf , or in other form 

—Au = f 

which is called the Poisson equation. If / = 0, it is called the Laplace 
equation in Q, which is a typical elliptic partial differential equation of 
second order. In dimension one, the Laplace equation becomes 

u = 0, 

So u = ax + b is a linear function of the single variable x. However, if 
the dimension n > 2, it is not easy to solve the Laplace equation on 
domains coupling with its boundary condition. 

In the following subsection, we try to solve some simple pde of first 
order. 
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2.2. Solve pde of first order. 

The principle for solving a first order PDE is to reduce it into ODE 
along characteristic curves. So one should know how to solve an ODE. 
We shall review some ODE theory when it is necessary. 

We study a first order partial differential equation. Our question is 
to solve the following equation (FDE): 

u t + au x = f(x,t). 

Assume that we have the initial condition u = u (x) at t — 0. 

We introduce the C-curve C(t) (called the characteristic curve) by 
solving dx/dt = a. Then we have 

C := C(t) :x = at + A. 

When t — 0, we have x(0) = A. On the curve C, we have 

du/dt = u t + au x = f(x, t) = f(at + A, t). 

Using integration along the curve C(t), we have 

u(x(t),t)=u(A,0)+ [ /(or + A,r)dr. 

Jo 

Given the point (x,t) on the curve C(t). Then we have 

x(t) = x 

and we can find that initial point on the curve is x(0) = A = x — at. So 
by the initial value condition, we have u(A, 0) = uo(A) = uo(x — at). 
Inserting A = x — at into the expression above, we find (FF): 

u(x, t) — u (x — at) + / f(x + a(r — t),T)dT. 

Jo 

In case when / = 0, we have 

u(x, t) = Uo(x — at). 

The method above is called characteristic curve method in the litera- 
ture. One may study examples in our textbook. 

We point out the following comparison lemma: Let u be a smooth 
solution to (FDE) with the initial data u(0) = u$. Let to be a smooth 
function satisfying 

w t + aw x > f(x,t) 

with the same initial data as u. Then we have w > u onH x [0, T]. 

This lemma is based on a very simple observation from the formula 
(FF) that if M = and / > 0, then u > for all t. In fact, let 

g(x,t) =w t + aw x - f(x,t). 
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Then we have g(x, t) > and 

w t + aw x = f(x,t) + g(x,t). 

Let v = w — u. Then we have 

v t + av x = g(x, t) > 

with v — at t — 0. Using the formula (FF) for v, we know that v > 0. 
Hence w > u. 

Usually, the comparison lemma is a very important property for 
PDE's. 
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2.3. A nonlinear PDE example. 

We have just used The characteristic curve method to study the 
linear pde of first order. We show here that this method can also be 
used to study the nonlinear pde of first order. By definition, a pde 
(PDE) 

P{u,Vu,D 2 u,...) = 

is called a linear equation if, for any given two solutions u and v of 
(PDE), au + v is also a solution of (PDE) for every constant a G R. 
Otherwise, it is called a nonlinear PDE. 

Let's study the following nonlinear PDE of first order: 

u t + (l + u)u x = - 

with the initial value condition u(x, 0) = \fx for x > 0. 
We first find the characteristic curve x(t) by solving 

tt =1 + 

Note that along the characteristic curve, we have 

du 1 
~dt ~ 2 

and 

u(x(t),t) = f -+ u(x(0),0) =*- + v^(0) 
On one hand, we set x(0) = A. Then we have 

u(x(t),t) = l - + VA 
Solving \fA from this relation, we get 

VA = u(x(t),t)- t -. 
On the other hand, we find the curve satisfying 

x(t) = A + t+ [ u(x(r),r)dT = A + t + / (- + \[A)dr. 
Jo Jo 2 

Then we have 

t 2 t 
x(t) = A + t + - + \fAt = t + (- + ^Af = t + u(x(t), t) 2 . 

We now solve u along the curve x(t) and get 

U (x(t),t) = y/x(t) - t. 
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Given any point (x,t) on the curve x(t), we have x(t) = x. Hence, the 
solution is 

u(x, t) = y/x — t. 
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3. LECTURE THREE 

3.1. Duhamel principle. 

Consider (ODE) 

u + p(t)u = 0, 

with initial condition u — f T at t — r. Then its solution to (ODE) is 
(S): 

u(t, t) = f T exp(- J p(s)ds). 

We now try to solve the equation (F): 

u +p(t)u = f(t) 

with initial condition u — at t — 0. Let f T = /(r) and solve (ODE) 
at t = t as above to get the solution (S). Then the function defined by 

u(t, r)dr 



is the solution to (F). This kind of method is called the Duhamel prin- 
ciple. 
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3.2. First order ODE and a comparison lemma. 

For the given ODE 

u' = a(t)u + b{t) 

we can write it as 

( e -/ to M u )' = e-/ to M&(*) 
So it has a general solution 

/t 
dse-f Sa{T) b(s). 

If «(0) = 0, we have (F): 

ft 



t{t) = e r^j dse -r. Msh 



One often use the argument above to prove the famous Gronwall 
inequality: Let (3 > and let u be a continuous function on [0,T] 
satisfying 

u(t) <a + 13 f u(r)dr. 
Jo 

Then we have u(t) < ae 131 on [0,T]. Hint: Let 

f[t) =a + /3 f u{r)dT. 
Jo 

Then we have /' = (3u. By assumption, u < f, we get 

f < Pf- 

Let (3{t) and f(s) be two smooth functions. We consider the following 
nonlinear ODE 

«' = /?(*)/(«) 
with initial data u = a at t = 0. 
Assume that 

«>' > /3(0/H 

with initial data w — a at t — 0. Then we have u>(£) > -u(t) for all t. 
Let 

g(t)=v/-0(t)f(w)- 

Then ^(i) > 0. Set 

v = w — u. 

Then we have 

v' == (3(t)(fW ~ /(«)) + »(0 = P(t)B(t)v(t) + 
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with v — at t — and 

B(t) = J f'(zw + (l- z)v)dz. 

Let a(t) = P(t)B(t) and b(t) = g(t) in (F). Then we have v > 0. That 
is w > u for all t. 

We remark that one can use Sturm-Liouville theory and variational 
principle to set up comparison lemma for second order ODE. 
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3.3. One dimensional wave equation. 

We shall follow our textbook [1] for the explanation. We consider 
the one-dimensional wave equation (OW): 

u u - u xx = 0, x e R, t > 0, 

with initial data 

u\ t=0 = u 

and 

u t \ t =o = Ui. 
Let's first assume that Uq — 0. Note that 

%-% = (d t -d x )(d t + d x ). 

Then we reduce the problem into two pde of first order. Let 

v = d t u + d x u. 

and 

dtv — d x v = 0. 

Solving u we get 

(x, t) = / v(x + (r — t), r)dr. 
Jo 

Using initial condition we have v(x, 0) = u\{x). We now solve v: 
v(x, t) = vo(x + t) — ui(x + t). 

So, 

u(x,t) = / u±(x + (r — t) + t)cLt. 
Jo 

Changing the variable 

s = x + 2r — t, 

we can get 



u(x,t) = -[/ ui(s)ds] := M U1 . 
If -Uq is non-trivial, but Mi = 0, then we can verify that 



u(x,t) = ^M U0 (x,t) 

is the solution to 

Utt ~Uxx = 

with initial data 

u\ t =o = U 

and 

Ut\t=o = 0. 
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This is the Duhamel principle for the wave equation. 

Hence, using the Duhamel principle, we can get the D'Alembert for- 
mula (EW): 

Y rx+t 

u(x,t) — -[(u (x + t) +uo(x — t) + / u 1 (s)ds]. 

* Jx-t 

for the solution to 

Utt ~Uxx = 

with the initial data 

u\ t=Q = u 

and 

u t \ t =o = Ui. 

From this formula we know that the solution to the problem (OW) is 
unique with the formula (EW). 

One can show that, for f s (x) = f(x, s), 

ft 



u(x,t) = / Mf s (x,t — s)ds 
Jo 



is the solution to 

u u - u xx = f(x,t) 

with the initial data 

u\t=o = 

and 

u t \t=o = 0. 
We leave this as exercise for readers. 

Remark One: We remark that if u = and u\ > 0, then u > for 
all t. For each t, \u(x,t)\ < M^+tMi is bounded provided |«o(^)| < M Q 
and < Mi. Similar bound is also true for u x and u t at each t 

provided |«ox(^)| < Co and |wi(x)| < C\. 

We now assume that 

Mo = 0, U\ = 

for all |x| > -R for some R > 0. Then from the expression (EW), we 
have 

u — 0, tx t = 

for |x| > -R + t. This is the finite propagation property of the wave 
equation. 

Remark Two: Define the energy for u by 



u\ + u 2 x )dx. 
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Then we have 

d f 

— E(t) = J (u t u tt + u x u xt )dx. 
Using the integration by part to the last term we get 

dl E ^ = J Ut ^ Utt ~ Uxx ^ dx = °" 
So the energy E(t) is a constant, which is 



2 

Remark Three: We now study the one dimensional wave equation 
on half line (HWE): 

u u - u xx = 0, x e R+,t > 0, 

with initial data 

u\ t =o = u 

and 

u t \ t =o = Ui. 

Note that Uo(t + x) is the special solution to (HWE). So we can 
use it to transform the problem into the case when uq — 0. Then 
we can use the odd extension method on the x-variable to extend the 
function u\ and the equation on the whole line. Then we can use the 
D'Alembert formula to solve (HWE). One may also use this method to 
treat one dimensional wave equation on an interval, however, we have a 
more beautiful method — Fourier method/Separation variable method 
to study this case. See the next section. 

Remark Four: We point out that we can use the D'Alembert for- 
mula to treat some special case of higher dimensional wave equation. 
Here is the example: Consider the following problem for 2-dimensional 
wave equation 

Uu U xx ^yy 

with the initial data 

u\ t=0 = 4>{x) 

and 

u t \ t= o = ij>(y). 

Then we split the problem int to two one-dimensional problems: One 
is 

u tt ~ u xx — 

with the initial data 

u\ t=0 = 4>{x) 
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and 

u t \t=o = 0. 

The solution for this one dimensional problem is 

u(x,t) = -[(<j>(x + t) + <j>(x-t)] 

The other is 

u tt ~~ u yy = 

with the initial data 

u\ t =o = 

and 

u t \ t= o = i/j(y). 

Its solution is 

y+t 



u(x,t) = -[ ij(s)ds}. 
z Jy-t 



Hence, the solution to our original 2-dimensional problem is 

I rv+t 
u(x, y- 1) = -[(<j>(x + t) + (f)(x-t)+ / ^{s)ds}. 
1 Jy-t 
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4. LECTURE FOUR 

4.1. The method of separation of variables 1. 

The separation variable method for boundary value problems is look- 
ing for a solution as a sum of special solutions of the form X(x)T(t). 
Let u(x,t) = X(x)T(t) be a solution to the wave equation 

u tt = a 2 u xx 

on [0, 1] x [0, oo) with the homogeneous boundary condition 

u (0,t) = = u(l,t). 
The boundary condition gives us (BC): 

X(0) =0 =X(1) 
Using the wave equation we have 

T"X(x) = a 2 TX", 

or 

a 2 T~ X' 

The left is a function of variable t and the right is a function of variable 
x. To make them be the same, they should be the same constant. Say 
—A. Then we have one equation for T: 

T" + a 2 AT = 

and the other equation for X: 

X" + AX = 

with the boundary condition (BC). The latter is the simple case of 
the Sturm-Liouville problem. The key step in the Separation Variable 
Method is to solve this problem at first. 

Case 1: When A < 0, the general solution for X is 

X(x) = de^ 1 + c 2 e~^ x 

Using (BC) we have 

ci + c 2 = 

and 

cie^ + c 2 e~^ = 0, 

which implies that c\ = c 2 = 0. 

Case 2: When A = 0, the general solution for X is 

X(x) = ci + c 2 x. 

Again (BC) implies that c\ = c 2 = 0. 
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Case 3: When A > 0, the general solution for X is 
X(x) = C\ cos(\^\x) + C2 sin(y/\x). 
The boundary condition (BC) gives us that c± — and 

c 2 sin(-\/A) = 0, 

and then 

\ = k 2 ir 2 , k = 1,2,... 

Hence 

X(a;) = -Xfc(x) = Cfc sin(/c7ra;) 

for fc = 1,2,.... 

Returning to the equation for T(t) with A = k 2 n 2 , we get 

Tfc(t) = a*; cos(A;7rat) + 6^ sin(/c7rat). 

So 

C/fc(x, t) = [a*; cos(knat) + b k sin(knat)} sin(A;7ra;) 

is a typical solution to the one-dimensional wave equation with bound- 
ary condition. 

We remark here that a easy way to see A > is using the integration 
relation: 



A \X\ 2 = IX 



-'|2 



We now look for a solution of the form 

+oo 

u(x, t) = ^^[dk cos(kTrat) + bk sm(kiiat)} sm(knx) 
k=i 

to the wave equation 

u tt = a u xx 

on [0, 1] x [0, oo) with the initial condition 

u = (f>(x), u t = ip(x), t = 

and the boundary condition 

u (0,t) = = u(l,t). 

The reason for us to use such an expansion is the Fourier series theory. 
For general boundary conditions, the expansion can be used , but we 
have to invoke the Sturm- Liouville theory (see Appendix). 

We need to use the initial condition to determine the constant a k 
and bk- So we have 



+oo 

(x) = 

fe=l 



a k sin(A;7ra;) 
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and 

+00 

i/)(x) = knabk sin(A;7r:r). 
fc=i 

Using the Fourier series expansion for <j)(x) and ip we have 



a k = 2 J 0(0 sin(kn£)d£ 

and 

= ^ y ip{Qsin(kir£;)d£;. 

In this way, we get at least formally a solution to our problem. This is 
a basic method to solve linear pde. 

If the boundary condition is not homogeneous, saying 

u(0,t) = a{t),; u(l,t) = /3(t), 

we have to make a transformation 

u(x, t) = v(x, t) + a(t) + x(p(t) - a(t)) 

and reduce the problem for v with the homogeneous boundary condi- 
tion. 



EQUATIONS FROM MATHEMATICAL PHYSICS 31 

4.2. The resonance problem. We now briefly discuss the resonance 

problem of the one dimensional wave equation. Consider 

u tt — a 2 u xx = Asm(ut), (x,t) € [0, L] x (0, oo) 
with homogeneous initial and boundary conditions 

u(0,t) = u(L,t) = 

and 

u(x, 0) = = u t (x, 0), 

where u ^ for j = 1,2,.... 
Then we have 

u{x,t) = y, — — sm(— — ) / sin(a;rj sm( — — [t — T))ar 
, a 3K L J L 

a i 



j 



EClj . t t . , . J 7TX . 
7rr{U Smud — UJj Smut) Sill — 



3 



where 
and 

Because we have 



a,- = — / Alrr) sin — — )<±c. 
LJ L 



2 2 



which may cause the coefficient 

a 



" oT — >■ oo. 



We say that we may have the resonance phenomenon caused by a;. 
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4.3. Wave equation in two dimensional disk. 

We denote B R (0) C R 2 the ball in the plane R 2 . Let's consider the 
2-dimensional wave equation 

U tt = a 2 (U XX + Uyy) 

on B R (0) x [0, oo) with the homogeneous boundary condition 

u(x,t) = 0, x e dB R {0) 
and the initial conditions 

u\ t =o = u 

and 

u t \ t =o = Ui. 

Introduce the polar coordinates (r, 9) by 

x = rcos9, y = r sin 9. 

Then 

Au = u rr + r~ l u r + r~ 2 U00. 

The reason that we use polar coordinates here is because of the bound- 
ary condition in our problem being on the circle. Let u = T(t)U(x,y) 
be the solution to our wave problem. Then 

T" + a 2 AT = 

and 

AU + \U = 0, U(R,t) = 0. 

It is easy to see that A = k 2 > 0. Setting U = h(6)B(r) in the latter 
equation, we have 

h"(6) r 2 B" + tB' + k 2 r 2 B 
h{9) ~ B ' ' 

Then we get two problems: 

h" +(3 2 h = 0, h(9 + 2n) = h(9), 

and 

B" + r~ 1 B > + (k 2 - f3 2 r~ 2 )B = 0, B(R) = 0. 
The last equation is called the Bessel equation. As before, we find that 

P = p n = n 

and 

h n (9) = A n cos(n9) + D n sin(n9). 
Inserting P — n into the Bessel equation, we have 

B" + r- x B' + (k 2 - n 2 r- 2 )B = 0, B(R) = 0. 



EQUATIONS FROM MATHEMATICAL PHYSICS 



33 



The general solutions of this equation are called the Bessel functions 
of order n. Just like trigonometric polynomials, Bessel functions form 
an orthonormal basis, and then we can use this basis to solve the wave 
problem above. One may look for reference book for more about the 
Bessel functions. Some people would like to define the Bessel functions 
of order n as the coefficients of the Fourier expansion: 

oo 

exp{irsm.6) = ^ B n (r)e me , 

n=— oo 

where % = \/— T and e md = cos(n9) + i sin(n0). Hence, 

1 f 27r 

B n ( r ) = — / exp(ir sin 9 — in9)d6. 
2tt Jo 

We remark that, if the domain is not a ball but a rectangle, we just 
let 

U(x,y)=X(x)Y(y). 
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4.4. * Wave equation in R 3 . 

The other way to understand the higher dimensional wave equation 
for us is to try the Duhamel principle to the wave equation in R™ 

u tt -Au = f; for (x, t) G R n x R + 

with the initial condition 

u = (f>(x), Ut — 0, at t = 

by assuming that we know the solution u(x, t) to the homogeneous 
equation 

u a - Am = 0; for (x, t) G R™ x R + 
with the initial condition 

u — 0, Ut — ip{x), at t = 0. 

We then use the descent method to study 

u tt = Am 

in R n x R in one dimension higher space R n+1 x R with initial conditions 

u — 0, Ut — ip(x), at t = 0. 
In dimension three, we first study the wave equation problem (3W): 
u tt - Am = 0; R 3 x {t > 0} 

with the initial data 

u — 0, -Uj = ip(x), at t = 0. 

We find that 

= ~a — 27 / i , {^)da x , 
A ^ a t Js at ( P ) 

where S at {p) = {x G R 3 ; r = \x — p| = at}. 
In fact, let 



I(x,r;g) = — / uda 

471_ JdB r (x) 

for g = g(y) G C 2 (R 3 ) and 

J := J(x, r; t) = rl(x, r, u) 

for u — u(y, t) be the solution to the problem (3W). Changing variables, 
we find that 



1 f 

I(x,r;g) = — / u(x + r6)d6 
^J{\e\=i} 
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and 

/ g(x + z)dz 



We now compute 



So, we have 



Note that 



and 



Jb r (o) 

= / r 2 dr / u(x + r9)d9 
Jo J{\e\=i} 

f R 

— — I 47rr 2 I(x,r;g)dr. 
Jo 

A x / Anr 2 I(x,r;g)dr 
Jo 

= A x g(x + z)dz 

Jb r (o) 

= / A x g(x + z)dz 
Jb r (o) 

= / A z g(x + z)dz 
Jb r (o) 

= / d R g(x + z)da z 

JdB R (0) 

= [ d R g(x + R9)R 2 d9 

JdBt(0) 

= =4nR 2 d R I{x,R;g). 

r R 

A x r 2 I(x,r;g)dr = R 2 d R I(x,R;g). 
Jo 

d R A x / r 2 I(x,r;g)dr 
Jo 

f R 

= A x d R I r 2 I(x,r;g)dr 
Jo 

= R 2 A x I(x,R;g) 

d R [R 2 d R I(x,R;g)] 
= 2Rd R I(x,R;g) + R 2 d 2 R I(x,R;g). 
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Hence, 

R 2 A x I(x, R; g) = 2Rd R I(x, R; g) + R 2 d 2 R I{x, R; g). 

Dividing both sides by R, we have 

A x RI(x,R;g) 
= = 2d R I(x,R;g) + Rd 2 R I(x,R;g) 
= d 2 R (RI(x,R;g)). 



Now, 



d 2 J R 
hi 

a 2 R 



dt 2 4tt .i dBl{0) 



/ utt(x + R9,t)d9 

JdBi (0) 



idl 



47r ./SBi(0) 



A x u(x + R0,t)d9 



a 2 r 

—A X R / u(x + R9,t)d9 

4?r JdBiiO) 



4tc~ 

= a 2 A x (RI(x,R;u)) 
= a 2 A x (RI(x,R;u)) 

„2' 



,8 2 J 



dR 2 ' 

This is a one dimensional wave equation. We now determine its bound- 
ary condition and initial condition. Clearly, we have the boundary 
condition 



At t = 0, 



and 



Then we have 



J(x,R;t)\ R=o = 0. 



J(x,R;t)\ t=0 = RI(x,R;0) = 0, 



J t (x,R;t)\ t=0 = RI(x,R;*P). 



■y rat+R 

J(x,R;t) = — / rl(x,r;t/j)dr. 

2o Jat-R 
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Hence, 



/ ,• J(x,R;t) 

u(x,t) = hm 



R^o R 

at+R 



= lim — — / rl(x,r:vb)dr 
R^2aRj at _ R K 

= -[RI(x,R;^)}\ R=at 

CL 

= tl(x,at;ifj) 
t 



ip((x + at6)d6. 
asi(o) 



That is, 



47T 



U 0M) = J—^ t I Hy)d°y 



'Sat(x) 

Then we use the Duhamel principle in the study of the wave equation 
■an -Au = f; R 3 x {t > 0} 

with the initial data 

u = (f>(x), u t = ip(x), at t = 0, 
and we have the famous Kirchhoff formula for the solution: 

u(p,t) = -^-r / r~ 1 f{x,t-a~ 1 r)dx 

47rfl JB at (p) 

Sat( P ) 

l ^a 2 t JSat{p) 
where 

B a t(p) = {x E R 3 ; r = |rr — p| < at}. 

From this formula, we can see the Huygens principle that the wave at 
x at time t — influences the solution only on the boundary of a cone 
originating from xq. We just remark that using the descent method, 
one can get a solution formula for two dimensional wave equation. 

One may see other textbooks for the other understanding of the 
Kirchhoff formula. 



+ j^t / i>{x)do- x 

^ aH Js at (p) 

+ ^7-^7 / <P{x)da x 



: J ,8 
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5. LECTURE FIVE 

5.1. The method of separation of variables 2. 

We begin with a typical example for heat equation with fixed bound- 
ary conditions. 

We now describe the main idea. Let u(x, t) = X(x)T(t) be a solution 
to the heat equation 

Ut = a 2 u xx 

on [0,/] x [0,oo). 

Then formally, we have 

a 2 T~ X' 

The left is a function of variable t and the right is a function of variable 
x. To make them be the same, they should be the same constant. Say 
—A. Then we have 

T' + a 2 AT = 

and (EX): 

X" + AX = 0. 

As in the wave equation case, we impose the homogeneous boundary 
condition 

X(0,t) =0 =X(l,t) 
to solve X and A. Multiplying both sides of (EX) by X and integrating 
over [0,1], we have 

A fx 2 = [ \X'\ 2 > 



So A > 0. Then we have 

X(x) = Xk(x) = Cksm(knx) 

for k = 1,2, ... and 

A = A fc = k 2 n 2 . 

Then we solve for T and get Tk(t) = T (0)exp(—a 2 X^t) by imposing 
initial conditions for T. The lesson we have learned is that the solution 
to heat equation has exponential decay to in t- variable as t — > oo. 

Using the superposition law, we can solve the heat equation by look- 
ing at solutions of the form 

A k exp{— a 2 k 2 n 2 t) sin(A;7r:r) 

k 

with the initial data 



4>{x) = A k sin(A;7rx) 
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Assume that we are given the heat equation 

u t = a 2 u xx 

with non- homogenous boundary condition, we need to use transforma- 
tions to make the boundary data into a homogenous condition, saying 

u(0,t) = = u x (L,t). 

Then we have to meet the following non-homogenous heat equation 

u t = a 2 u xx + f(x,t). 

To study this problem, we need to solve for {Xj} and A = Xj to the 
Sturm-Liouville problem 

X" + AX = 0, 
with the boundary condition 

X(0) = = X x (L). 
Once this done, we do the expansion for / 

and look for solution u(x,t) such that 

3 

Then we reduce the non-homogenous problem into the ODE 

u. = - a 2 \ jUj + fj(t) 
with suitable initial condition. Hence, we have the expression 

Uj (t) = e- A ^( % (0) + fe^f^dr 

Jo 



and 

u(x 



,t) = 5>-^K(0) + fe^MrWXjix). 

Jo 



Here Uj(0) are determined by 



We omit the detail here. 
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5.2. Energy bound. 

Let L be a partial differential operator on the domain D C ~R n such 
that for some > X 3 ;6 R, we have 

with a fixed homogeneous boundary condition. Define the function 
space E such that the functions {<f)j} form an orthonormal basis for 
the space E with its norm given by 

for u = a,j(j)j(x) G E. 

Given a smooth function / e C 1,2 {D x [0, oo)). Assume that the 
expansion u = Uj(t)(f)j(x) is a solution to the partial differential 
equation: 

u t = Lu + /, t > 0, 

with the initial data u — at t — 0. Hence, %(()) = for all j. 
Assume that the known function / has the expansion 

with its norm on D x [0, oo) defined by 

/oo 
|/,-(f)| 2 df<oo. 

3 - 

Then we have 

u J (t) = \u ] (t) + f)(t). 
Hence, as before, we have the expression 

Jo 

Note that 

K(0 - V,-(t)) 2 = «i(*) 2 + \-«i(*) 2 - 2X j u j (t)u' j (t) = /,(*) 2 . 
Then we have 

Tkw 2 + A; Ui (o 2 )dt - a^-(t) 2 = / T /,(t) 2 . 

Jo Jo 
Note that we have assumed that Xj < 0. Hence, we have 

/■oo poo 

/ («;(o 2 +A;«,-(o 2 )cft< / /,w 2 , 

Jo Jo 
which implies that 

IM 2 + A 2 N|<||/|| 2 , 
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where A 2 = infj{\^}. This is the energy bound for the solution as we 
wanted. 
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5.3. Fundamental solution to Heat equation. 

We study the heat equation (HE) on R n x [0, +oo): 

u t = Au, 

with the initial data u = g(x) at t — 0. 

Assume that u is a smooth solution of (HE). Note that for A > 0, 
the function 

u\(x, t) = u(Xx, \ 2 t) 

is also a solution of (HE). 

So we try to look for solutions of the form 

u(x,t) = tPf(\x\ 2 /t). 
Set r = \x\ 2 /t. Then we find that 

d 
d\x\ 
= ^xlt^f 

d 2 
o\x\ z 

= 2^- 1 /' + 4|x| 2 ^- 2 /" 
= 2t p ~ l f + 4rt /3 - 1 f" 



n — 1 / 

Au = u H — — u 

\x\ 



So, 



2^-y+ 47^-7" +2(n- 

4r^- 1 /" + 2n^- 1 /'. 



-^/ + /'+4/" + -/' = 0. 
r r 



Let P = —n/2. Then we have 

4(r n/2 /')' + (r n/2 /)' = 0. 
Then for some constant A, 

Ar n/2 f + r n/2 / = A. 
Take A = 0. Then we have 

4/' = -/. 
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Assume that /' — > C — const, as r — > +00. Then we have / — > as 
r — > +00. So 



Choose 5 = l/(47r) n / 2 . Define, for i > 0, 



1 



This is called the Gaussian kernel, which is also called the Poisson 
kernel. It is easy to verify that 



Then u(x,t) satisfies the heat equation (HE) and as t — > 0, u(x,t) 



Roughly speaking, the theorem says that this u(x,t) is the solution 
to our problem at the beginning of this subsection. 

Proof. For t > 0, we can take derivatives into the integrand term, and 
then u satisfies (HE) since K does. 

Without loss of generality, we now assume n — 1. Assume that there 
is a positive constant M > such that 



We use two step argument. 

Step one: Given an arbitrary small e > 0. For any fixed point x G R, 
we have some So > such that, for \x — y\ < So, 



Step two: We take Si > small enough such that for < t < Si, 




Theorem 2. Given g(x), we let 




g(x)\ < M, x e R. 



g(y)-g(x)\ < -. 
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Then we have 



\u(x,t) -g(x)\ = | J K(x -y,t)(g(y) - g(x))dy\ 

< I / K(x-y,t)(g(y)-g(x))dy\ 

J{\x-y\>S } 

+ I / K(x-y,t)(g(y)-g(x))dy\ 

J{\x-y\<S } 

< 2M K(x-y,t)dy 

h\x-y\>to} 

+ e K(x-y,t)dy 

JUx-v\<Sn} 



'{\x-y\<S } 
< 2M ' , • ' 6. 

AM 2 



Here we have used the fact that K > and 

J K(x-y,t)dy = l 
for any t > 0. This says that u(x, t) — > g(x) as t — > 0. 



□ 



We now study the behavior of the solution defined above as \x\ — > oo. 
Theorem 3. Assume g G C°(R n ) satisfies 

g(x) — >■ 0, as |x| — > oo. 

T/ien for the solution u(x,t) defined above, we have 

u(x,t) — > 0, V x m(o;, t) — ► 0, as \x\ — > oo. 

/or eac/i t > 0. 

Proof. Fix £ > 0. Choose i? > sufficient large. Then by our assump- 
tion on g(x), we have 

s«PB H (x)k(aj)| -»■ 0, \x\ -> oo. 
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Note that 



\u(x,t)\ < / K(x - y,t)g(y)dy 

< sup BR (x)\g(x)\ / K(x-y,t)dy 

JB r (x) 

+ sup R n_ BR{x) \g(x)\ / K(x-y,t)dy 

JR n -B R (x) 



poo 

< sup BR{x) \g(x)\+C(t)sup R n\g(x)\ / e~ p2/At p n 

JR 

-> + o(R) 

as \x\ — > oo. Here C(t) is a constant depending only on t and 
as small as we like for R large. Hence, we have showed that 

\u(x, t)\ — > 0, 

as \x\ — > oo. 

Using the polar coordinates (p, 0), we can also show that 

\Vju(x,t)\ -»• 0, 
as b| — > oo. We omit the detail here. 
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5.4. Comments on Fundamental solutions. 

The Gaussian kernel K(x — y, t) satisfies the heat equation for t > 0. 
For t — 0, we understand it in the generalized sense that 



That is, we consider 5 x (y) as a function defined on the function space 
C^{R n ). Then, S x is a linear functional on C °°(R). This 5 X {-) is called 
Dirac's delta function or Dirac measure. One can also consider the 
delta function 5 (x) as the limit of a sequence of nice functions. For 
example, assume that n = 1 and let 



Jn 

and for s^O, ^(ic) — > and at x = 0, ?j(0) — > oo; or for e Cq°(R), 



as j — > oo. If so, we define /img^ = <5 . We ask the reader use the first 
concept of limit to understand delta function. I would like to point 
out that the latter limit is more rigorous definition. We may find other 
sequences which have 5q as its limit. For example, let 



p(x) = me i-M 2 , for \x\ < 1, p(x) = 0, for \x\ > 1 
where m > is a constant such that 



K(x-y,0) = 5 x (y)- 
Here we define S x (y) as the generalized function by 

< Sx(y),g(y) >■= S x {y)(g(y)) = g(x). 



Then 





= < S , > . 




and let 



Pt(x) 



1 



p(x/t), 



t 



where t > 0. Then 




and Zim 
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Generally speaking, for any given continuous function u(x) in R n , 
we can also consider u as a linear functional on the space C^°(R n ) with 

<u,(f)>= / u(x)<fr(x)dx. 
Jr.™ 

With this understanding, we say K(x — y, t) is the fundamental solution 
the heat equation on R n x R + . We emphasize that we call K(x — y, t) is 
a fundamental solution to the heat equation on the whole space because 
K(x — y,t) satisfies the heat equation fort > 0, and we consider K(x — 
y, t) as a generalized function which is the delta function for t — 0. 

We now consider the heat equation on the half line [0, oo) with the 
boundary condition 

u(0,t) = 0, 

and the initial condition 

u = <f>(x), t = 0. 

For this case, the fundamental solution is 

T(x, y, t) = K(x -y,t)- K(x + y,t). 

So in some sense, we use the odd extension to keep the boundary 
condition. So the function defined by 

u(x,t) = / T(x,y,t)g(y)dy 

satisfies the heat equation on the half line [0, oo) with boundary con- 
dition u(0, t) = and the initial value condition 

u = <f>(x), t = 0. 

If we consider the heat heat equation on half line [0, oo) with the 
boundary data 

u x (0,t) = 0. 

Then the fundamental solution to his problem is 

F(x, y, t) = K(x -y,t) + K(x + y,t). 
Here we have used the even extension to keep the boundary condition. 
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6. LECTURE SIX 

6.1. Fundamental solution to Laplace equation. 

Consider the Laplace equation (L) on R n : 

Au = 0. 

A solution to the Laplace equation in a domain is called a harmonic 
function. 

Let r > such that 

i 

Then we have 

rdr = x l dx l . 

i 

Then 

dr x % 
dx % r 

Let u = v(r) be a solution of (L). Note that 



and 

Then we have 
Solve 

We find that 



r 

\2 



,.1 (x 1 ) 2 . „{x l ) 
U x i x i =v( — ) + V 



/y /y 3 /y 2 



// n — 1 / 
Au = v H v . 

r 

i, n — 1 / 
v H v = 0. 

r 

/ c 
v = 



j.n—1 



for some constant c and 

f = alogr + 6, /or n = 2, 

and 

f = ar 2_n + 6, /or n > 3, 
for some constants a and 6. 
Choose 6 = 0. Define 



and 



T(x) = — logr, for n = 2, 

27T 



r(x) = — ^ lry ^ r 2 ~ n , for n>3. 

y ' n(n-2)\S n - 1 \ J 
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Here IS" 1-1 ) is the volume of the sphere S 71-1 . 

Let T(x,y) = T(x — y) for x,y e R n . This function T(x,y) is called 
the fundamental solution (or Green function) to the problem (L). In 
fact, we can consider 

-AyF(x-y) = 5 x {y). 

If so, we have 

u(x) = <5 x (y),u(y)> 

= - < AT(x-y),u(y) > 

= - AT(x - y)u(y)dy 
V{x - y)Au{y)dy 



L 



R™ 

F(x -y)f(y)dy. 

/R" 

To understand this, we need to use the second Green formula. Here we 
recall the second Green formula as follows. Given two smooth functions 
defined on a bounded domain D with piece-wise smooth boundary. We 
denote v the outer unit normal to the boundary dD. Then we have 

/ (vAu — uAv)dx = I (vVu ■ v — uVv ■ v)do~. 

JD J 3D 

Theorem 4. Given f e Co(R"), i-e, f is continuous in R™ and f(x) = 
in the outside of some large ball. Let 



u(x) = T(x- y)f(y)dy. 

Then the function u satisfies the Poisson equation 

-Am = /, in R n . 

Proof. For simplicity, we assume n — 2. Take R > 1 large such that 
/ = outside Br. Fix x G R 2 . Using the second Green formula on 
the region B R — B € (x), where e — > 0, we get that 



T(x,y)(-Au(y)) 
lim / F(x,y)(-Au(y)) 
lim / (r(x,y)V r u(y) 

e ^JUv-x\=e} 



'{\y-x\=e} 

V r r(x,y)f(y))da, 
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where r = \x — y and V r u{y) = §^(y)- 



Note that 



/ 



\T(x,y)V r u(y)\do 



l{\y-x\=e} 
1 f 1 

< — / ln-\V r u{y)\ 

<{\y-x\=e} r 



2ir 



< sup\ Vu\e\lne\ 
- 0, 



and 



1 

2n 



V r T(x,y)u(y)da 

{\y-x\=e} 



V r ln(r)u(y)da 

{\y-x\=e} 



= ^ / -f{v)d° 

27T J{\y-x\=e} r 

= — / «(y)d<T 



— > w(a;) 

as e — > 0. This shows that the function u satisfies 
-/ F(x,y)Au(y)dy = u(x). 

J-R? 

We ask our readers to verify the n > 3 case. □ 

Let's study the Poisson equation on bounded domain. Let D C 
R™ be a bounded domain with piece-wise smooth boundary. We now 
consider the problem (PD): 

—Am = /, in D 

with the boundary condition u(x) = <f>(x) for x G dD. Then using 
the second Green formula and the fact v = T being the fundamental 
solution on the whole R™, we get (Ul): 

u(x) = / T(x-y)f(y)dy 
Jd 

+ [ (F(x - y)Vu{y) ■ v - 0(y)Vr(x - y) ■ v)da y . 

J 3D 

This expression still contains a unknown term for Vu(y) ■ v on the 
boundary dD. To overcome this drawback, we introduce, for each 
fixed point x G D, the harmonic function g(x, •) with its boundary 
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data g(x,y) = —T(x — y) for y G dD. Using the expression above for 
v = g(x,y), we have (U2): 

= / g(x,y)f(y)dy 
Jd 

+ / (g(x, y)Vu(y) ■ v - <Ky)Vg(x, y) ■ v)da y . 

JdD 

Set 

G(x, y) = T(x -y)+ g{x, y). 
Note that G(x, y) has zero boundary value for y, i.e, 

G(x,y) = 0, yedD. 

Adding (Ul) and (U2), we obtain that 



(x) = [ G(x,y)f(y)dy 
Jd 

- / <t>{y)VG(x,y) ■ vda y . 

JdD 



We call G(x, y) the Green function (or fundamental solution )for the 
problem (PD). 

Let's consider the boundary value problem for Laplace equation. 
Assume that u G C 2 (D) n C(D) satisfies 

-Am = 0, in D C R™, 

with the Dirichlet boundary condition 

u{x) = <j>{x) e C(dD), x G dD. 

As before, for x G D, we have 

u(x) = / F(x - y)(-Au(y))dy 
Jd 

[ (F(x - y)Vu{y) ■ v - (f)(y)VT(x - y) ■ v)da y 

JdD 



+ 



idl 



(p(y)VT(x - y) ■ vda y . 

IdD 

Let 

P(x,y)=VT(x-y)-is 

for x G D and y G dD. We call P(x,y) the Poisson kernel to the 
boundary value problem in D. Then it can be verified that the function 
defined by 



u\x) = 



P(x,y)<f>(y)do-y, 

dD 
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which is called the Poisson formula, is the solution to the boundary 
value problem in D. 

We now point out a nice observation about the construction of har- 
monic functions. Given two harmonic functions U\ and u 2 on the same 
domain D C R n . We want to find a harmonic function U(x,y) defined 
on the one dimensional higher cylindrical domain D x [L 1 , L 2 ] such that 




U(x, L 2 ) = u 2 (x). 
This function can be achieved by setting 



U(x,y) = ^ — y-M x ) + 



V-Li 
L 2 — L\ 



u 2 (x). 
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6.2. Mean value theorem for harmonic functions. 

Given a harmonic function u in a domain D. We take a ball Br{j>) 
in the domain D. 

Note that x = p + rd, where r = |x — p\ and = Using the 

divergence theorem, we get 



means the radial derivative of function u taking value at x = p + tO. It 
is not a composite derivative for independent variables (r, 9. Hence, 




Then we have the the right side of the identity above is 




Here, we want to point out that the notation 



du . 




This implies that 




is a constant for all t < R. Then 




Using the mean value theorem and sending t — > 0, we have 




Then 




If we take integration in R, we have 




Hence, we have 
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Theorem 5. Given a harmonic function u in a domain D. We take 
a ball Br{p) in the domain D. Then we have 



and 



1 f 

U(p) = r— . . . / Udd. 

yp> \dB R { P )\ J 9Br(p) 



u(p) = / udx. 



\Br(p)\ 



Br( P ) 



Mean value theorem is a basic property for harmonic functions. As 
an application we prove the following result. 

Theorem 6. (Liouville Theorem) . Ifu is a bounded harmonic func- 
tion in the whole space R™. Then u is a constant. 

Proof. Method One. Let M > such that \u(x)\ < M for all x E R n . 
Note that 

1 



u ( x ) = i R 77vH / u(x + y)dy. 
\Br{V)\ J Br (o) 



Hence, 

d xi u(x) = 1 / Viu(x + y)dy 
\Br{v)\ J Br (o) 



1 



\B R yv)\ j 9Br (o) 
where v i = x l /R, and then, 

M\dB R 

\d x iu(x)\ < 



b r (o) 

h>iu(x + y)d<jy, 



\Br(0)\ 

as R — > oo. This implies that tt is a constant. 

Method Two. Fix any two point x and Xi. Then for any R > 0, we 
have 



u(x ) - = 1 ( / udx - / ttob) 

= [RTrnT ^ / udx- I udx). 

I^M U JI ./Br(xo)-.B,i(xi) J B R (xi)-B R (xo) 



Hence, 

|u(xo) — u(xi) < ( / |w|dx + / |u|dx) 



< 



2M| J B R (x 1 )- J B R (x )| 



15 
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as R -> oo (since \B R ( Xl ) - B R (x )\ « iT -1 and |5 R (0)| « iT). That 
is to say, u is a constant. □ 
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6.3. * Monotonicity formula for Laplace equation. 

This subsection is only written for readers not just for fun, but also 
for its fundamental role in the regularity theory of nonlinear elliptic and 
parabolic equations. It is not required for the examination for students 
even though it is really important for modern PDE. Recently I and 
my students have found a monotonicity formula for solutions to a very 
large class of elliptic and parabolic equations. We can not exhibit it 
here, however, we would like to show the classical monotonicity formula 
for harmonic functions. 

Assume that u is a harmonic function in a ball Br := Br{j>) with 
center at the point p and of radius R. It is easy to see that the integral 

/ \Vu\ 2 dx 
Jb r 

is a non-decreasing function of R. To our surprise, so is 

R 2 - n [ \Vu\ 2 dx. 
Jb r 

Why? let us prove it. Recall that 

Au = 

in the ball B R . Let X{x) = (Xj(x)) = x — p be a radial type vector 
field in Br. Recall that on 8Br, the outer unit normal is 

v = (x — p)/R- 

We denote by 

du _ d 2 u 
Ui =dx^ = V * M ' U - = dx^' 

We now have 

= -/ AuVu-X = - Vu-uVu-Xda+ I Vu- V{Vu-X)dx. 

J Br J 8Br J Br 

Then we re-write it as 

/ Vu ■ uVu ■ Xda = Vu- V(Vu ■ X)dx. 
JdB R J Br 

Note that 

/ Vu ■ vVu ■ Xda = R I \d r u\ 2 . 

JdB R JdB R 

and 

/ Vu-V(Vu-X)dx = / ViuViiVjuX^dx = / (\Vu\ 2 +u ij u i X j ). 

J Br J Br J Br 
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Using integration by part to the last term, we have 

J b UiiUiXj = l -J V j \Vu\ 2 X j = ^J gB |V*f ~\j n\Vu\ 2 . 
Then we have 



Vu • V(V« • X)dx = ^ / |Vm| 2 + - / |Vw| 2 , 



2 



9B R 



and 



(n-2) / 
Jb r 



V< = R I (\Vu\ 2 - 2\d r u\ 2 )da. 

dB R 



We now compute 

d R [R 2 - n I \Vu\ 2 } = R 2 ~ n ( \Vu\ 2 + {2-n)R 1 - n I |V 

J B R J 8B R J B R 

= 2R 2 ' n [ \d r u\ 2 



dB R 

> o. 

This implies that the function 

R 2 ~ n I IV// 1 



-/ 

Jb, 



'b r 

is a monotone non-decreasing function of R. More precisely, we have 
proved the following result. 

Theorem 7. Assume thatu is a harmonic function in a ball Br. Then 



R 2 



n ( \Vu\ 2 -r 2 ~ n ( \Vu\ 2 = 2( \x-p\ 2 \d r u\ 

J B ft J B r *J B R — Br 



' B R J B r J B R —B r 

for all r < R. In particular, the function 



r 2 ~ n I |Vd 



is a monotone non- decreasing function of r . 

Monotonicity formulas play a very important role in the regular- 
ity theory of weak solutions to nonlinear parabolic and elliptic partial 
differential equations/ systems. 



L( 
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7. LECTURE SEVEN 

From the previous sections, we know that solving a PDE is not easy 
matter. To handle PDE on unbounded domains, in particular, PDE on 
whole space, people try to use integration type of transformations to 
reduce ODE into algebraic equations and to reduce linear PDE into to 
ODE. The famous Laplace and Fourier transforms are for this purpose. 

In this section, we always assume that our functions have a good 
decay at infinity such that the integrals converge. 

7.1. Laplace transforms. 

We define, for a function y = y(t) G C°[0, oo), the Laplace transform 

by 

poo 

<(y)( 8 ) = / y(t)e- st dt. 
Jo 

Remark: 

It is clear that if there is some constant M > such that 

\y(x)\ < e Mx , forx>0 

then the Laplace transform for y is well defined. 
From the definition, we have 

MVi + a V2) = Hvi) + aL(y 2 ) 
for every real number a G R. 
Note that 

PCX) 

Uy')(s)= / y'(t)e~ st dt = -y(0) + sh(y)(s). 
Jo 

By induction, we get 

h(y^)(s) = s k h(y)(s) - Xy^oy-''. 

3=1 

We can use this property to solve the following ODE: 

y" + ay' + by = g(t), 

with initial conditions 

y(0)=0 = y\0), 

where a, b are constants and g(t) G C°[0, oo). In fact, let Y(s) = 
L(y)(s). Then we have 

s 2 Y + asY + bY = h(g). 

Hence, we have 

Y _ Hg) 

s 2 + as + b 
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Denote by Lr 1 the inverse of the Laplace transform. We remark that, 
using complex analysis and inverse Fourier transform (see next section) , 
it can be proved that 



where a is the real part of s = a + ib, which is a complex variable in 
the plane. Then we have 



In general, we have a list of Laplace transforms of special functions to 
use. So we shall invoke to special book on Laplace transforms. 
Example. We now compute one example to close this part. Given 




y(t) = t n . 



Then we have 



Jo 

We do the computation step by step. 





s n+l ■ 
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7.2. Fourier transforms. 

We first introduce the Schwartz space S. We say / is a Schwartz 
function on R if / is smooth and for any non-negative integers m, n, 

\x\ m \f (n) {x)\ ^0 

as x — > oo. Then the Schwartz space is the space of Schwartz function 
on R. Note that C£°(R) is a subspace of S. The important example is 
that the function 

g A (x) = e~ Ax2 

is in the space of Schwartz space, where A > is a constant. We denote 
g(x) = g±(x). 

Define the Fourier transform of a Schwartz function. Let / be a 
smooth function in S. Define 



/(0 = -±= J f(x)e-^alx, 



which is called the Fourier transform (in short, F-transform) of / in 
R. Set F(f) f. 

Remark: From the definition formula, we see that the Fourier trans- 
form makes sense for any function / G L 1 (R n ), that is, 

\f(x)\dx < oo. 



R" 

From the definition of the F-transform, we can see that the Schwartz 
space is the invariant set of F-transform, that is, 

F : S -> S. 

Example 1: Let's compute the Fourier transform of the function 
Qa^x) above. By definition, 



Note that 



g2(0 = -±= J < A ' , " v/.r. 



-Ax 2 - ix£ = -A(x + —Cf 



2A S ' 4A 2 ' 
and 



/ 



e Ax2 dx 



Then we have 
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From this we can see that, for A 3 — |, g& is the eigenfunction of the 
F-transform. 

For t > 0, let A = ^ and let = ^7^e™3^. Then 

F[^,0](0 = e-^ 2t . 

We shall use this expression. 
Example 2: Let 

qj(x) = — , for \x\ < and qj(x) = 0, for \x\ > 1/j. 

Then 

1 r 1 "-^ 1 sin « 



1 

2V27T.7 7-i/ 7 - 

Hence, 



£(0) = 



as j — > oo, and 



/2tt 

Given /i, ^ A G R. Let 

(r h f)(x) = f(x-h), (5 x f)(x) = f(Xx). 
Then by definition, we have 

(S7)(0 = e- i(fc€) /(0 

and 

(S7)(0 = A-7(e/A). 

Hence, for 

(M(^)=e- AV , 

we have 

Given two functions / and g. Define their convolution by 



f*g{x)= J f(x -y)g{y)dy. 
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Then we have 

ix 



F(f*g(x)) = -±=j f*g(x)t 

= = J e~ ix ^dx J f{x - y)g{y)dy 

= -J= J g(y)e-^dy J f(x - y^^dx 



'2tt/ • g. 

The most important property is that the F-transform of the derivative 
f x is a multiplication: 

LAO iCfii). 

This property gives us a chance to define the inverse operator d~ l of 
d x by the definition 

d^fiO = (iO-'ttO- 

Clearly, we have 

d-'dj = f. 

Remark: One may see more interesting properties about Fourier 
transforms in the subject called "Harmonic/Fourier Analysis" , which 
is one of the beautiful theory of modern mathematics. 

Let's going back. By an elementary computation, we have the fol- 
lowing symmetry formula: 



J u{x)g{x)dx = J u(y)g(y)dy. 



Define 

/(0 = -L / f(x)e^dx, 



which is called the inverse Fourier transform of /. In short we write 
it by F*(f). 

By an elementary computation and the symmetry property, we have 
Proposition 8. 

/ = / 

and 

/ = /■ 

Proof. Let's verify the latter equality. In fact, note that 

g A (x) -> 1, 
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as A — > and 



as A — > 0. Then, by symmetry property, 

/(0) = wJ nm 



\=lim A ^ J gl{y)f{y)dy 
= /(0). 

In general, we let fi(y) = f(y + x). Then 

fix) = aw = y A(0de = ^ / / v ^- 



Using this we can prove 
Theorem 9. (Plancherel Theorem): Given u G S 7 we have 

\u\ 2 dx = I \u\ 2 dt; 



Proof. In fact, using the inverse Fourier transform, we have that 
u(x) = --j= J u(y)e~ txy dy = u(x). 



Let g = u in the above equation. So 

u = g. 

Using the symmetry formula. Then we have 



\u\ 2 = uu = ug = ug = / \u\ 2 . 



One can define n-dimensional Fourier transform of / in R™ by 
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Here x ■ £ is the inner product in R n . The inverse Fourier transform of 
/ in R n is defined by 




Similar properties in dimension n are also true for n-dimensional 
Fourier transforms. In particular, Plancherel Theorem in R™ is true. 
One important fact lying behind the Fourier transform is that for the 
linear operator 

A[u] = u (x), 
we have the spectrum property: 

A[e tx Z] = -|£|V^. 
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7.3. Heisenberg uncertainty principle. 

For any smooth function u (maybe complex-valued) vanishing at in- 
finity with 

|2 



/ 



U\ 



we have 



(jx 2 \u(x)\ 2 dx)(je\u\ 2 dO>\, 

with its equality true only at 

u{x) = Ae~ lix \ 

where (3 > and A 4 = 2(3 /n. 

The inequality above is the so called Heisenberg uncertainty princi- 
ple. 

Proof. Integrating by part, we have 

l = y \u\ 2 = — 2 J xu ■ u < 2 J \x\\u\\u\. 
Using the Cauchy-Schwartz inequality we have (H): 

1<4(|xV)(|| M '| 2 ), 

with equality at 

u = Xxu. 
The latter condition implies that 



u 



Ae Xx2 ' 2 . 



Using J u 2 = 1 we have A = —2(3 and the relation between A and (3. 
Note that 



Inserting this into (H), we then get what wanted. 
Notice that, as a byproduct, we have 

l -ju 2 <{jx 2 u 2 Y' 2 {j\u\ 2 yi 2 , 

for every u G Cq(R). 

There are a lot of applications of Fourier transforms in the various 
branches in sciences and technologies. I hope our readers can read more 
books on Fourier transforms and Fourier analysis. Please read as more 
as you can. 



66 LI MA 

8. LECTURE EIGHT 

8.1. Application of Fourier transforms. 

First, let us use the F-transform to solve the wave equation 

u tt = Au, on R n x R + 
Take F-transform for the variable x on both sides. Then we have 

Utt(U) = -\tfu(U). 
Solving this ODE, we get 

^ = ^,0)cos(|e|t)+^(e,0)^p. 

Using the inverse F-transform we find the solution expression: 

u(x, t) = F*(u(Z, 0) cosflflf) + 0)^|p). 

Set 

u(x,0) = f(x), u t (x,0) = g(x). 
The solution can be written as 

Using the relations 



and 



we have 



cos(|e|t) = ^(e^ + e-^) 
sin(|e|t) 1 imt mt 

lei im [ 1 



where 



and 



u(x,t) = J (f + (0e^ x+m + /-(Oe i(€x - |€|t) )de, 

2 



Using the similar approach, we can use properties of F-transform to 
solve the heat equation 

u t = Au, on R™ x R + 
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Take F-transform for the variable x on both sides. Then we have 

This is a first order ode, and we can solve it to get 

u = u(^0)exp(-\£\ 2 t). 
Using the inverse F-transform we have 

u(x,t) = F*(u(Z,0)exp(-\Z\ 2 t)). 

Let 

g t (x) :=g{x ,t) = ^e-M a <« 

Then 

9t(Z)=exp(-\Z\ 2 t). 
Given the initial data u(x, 0) = f(x) in R n . Let 

Let 

H t (z):=H(z,t) = j^e-W a t« 

Then 

u(x,t) = f*H t (x). 

Hence, u{£,,t) = fexp(—\^\ 2 t), and Ut(£,t) = — |£| 2 w(£, t). This shows 
that u — f * H t is the solution to the heat equation with initial data 
u(x,0) = f(x). 
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8.2. Laplace equation in half space. 

Our problem is to solve the Laplace equation 

Am = u xx + u yy = 
on the half space = {(x, y); y > 0} with boundary condition 

u(x,0) = f(x). 

To use the F-transform method, we think of u y {x, 0) = 0. Then we 
have 

-\t,\ 2 u(t,,y) + Uyy(t„y) = o. 

Solve this and we find 

u(Z,y) =u(Z,fyexp(-\Z\y) = f(O ex p(-\Z\y)- 

Define 

P y (x)=F*(exp(-\-\y))(x). 
Then by an elementary computation, we have 

P y (x) = F*(exp(-\-\y))(x) 
1 i e-Wv+i**^ 



= — / 

V2^ Jo 

L_ [° 



This implies that 



Hence 



and 



1.1 1 . 
-1=( + -)■ 

V 27r y + ix y -ix 
1 2y 



P y (x) 



2ir x 2 + y 2 ' 



V ^7T 



«(z>2/) = -^=f* p v( x ) 

V 27T 



is the solution wanted for the problem above. 
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8.3. Derivative bound. We ask the following question. Given a 
smooth function u G S(R n ). Assume the bound M of \u\ and the 
bound M2 of |Au| on R n . Can we give a bound of |Vw| on R n in term 
of M and M 2 ? We answer this question below. 
Let 

-Am + u = f, in R n . 
Using the Fourier transform, we get 

Then we have 

Define 

B2(x) = J- F*(—^ )(x). 

Then we have 

u(x) = B 2 *f(x) 

and 

\Wu(x)\ = \VB 2 *f(x)\ < \VB 2 \*\f\(x). 
Using the knowledge from calculus we can show that 



/ \VB 2 \(y)dy:=B <oo. 



Notice that 

< max\u\ + max\Au\, 
where max means taking maximum value over the domain R n . Hence, 
we have 

|Vw(a;)| < B(max\u\ +max\Au\). 
Using the function u(rx) (with r > to be determined) to replace u 
in the inequality above, we have 

r\Vu(rx)\ < B(max\u\ + r 2 max\Au\). 

Since r > can be any positive number, we get 

\Vu(rx)\ 2 < AB 2 max\u\max\ Au\. 

Therefore, 

max\Vu\ < IB \/ max\u\max\ Au\. 
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9. LECTURE NINE 

9.1. Maximum principle (MP) for Laplace operator. 

The Maximum Principle (MP) is a beautiful property for solutions 
to Laplace equations and heat equations. One can use the Maximum 
principle to get uniform bound for a solution and its derivatives. 

Consider a smooth solution u to the following elliptic partial differ- 
ential inequality in a bounded domain D in the plane R 2 (IPDE): 

^xx ~t~ ^yy ^ 0* 

Note that there is no maximum point of u in the interior of the 
domain D. For otherwise, assume (xo,yo) is and then we have 

Vu(x ,y ) = 0, D T D T u(x ,y ) < 0, 

for any direction T e R 2 . The latter condition implies that the matrix 
D 2 u(x ,y ) is non-positive, that is, 

D 2 u(x ,y ) < 0. 

Hence, 

(u xx + u yy )(xo,yo) < 0, 
which contradicts to (IPDE). Therefore, we have 

supi>u = SWPqdU. 

Assume that we only have 

^xx ~t~ ^yy — ^ 

in the domain D. Then we choose 

w = u + ee 13 * 

for small positive constant e and large positive constant (5. Note that 

Ae px > 0. 

Then we have 

Aw > 0. 

Hence, we have 

supow = SUPQDW. 

Sending e — >, we get that 

supnu = supodu. 

Clearly, the above argument is true for general domains in R n . This 
is the maximum principle for elliptic pde of second order. 
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Theorem 10. Maximum Principle. Assume u is a smooth function 
satisfying 

Au > 0, in D. 

Then we have 

sup D u = supg D u. 

Remark: (1). It is clear from our argument that the function u G 
C 2 (D) n C(D) satisfying 

u xx u yy < 0, in D, 

can not reaches its maximum in the interior of the domain D. In fact, 
if it has the maximum point at (xo,2/o) i n the interior of D, then we 
have 

^xx — 0> ^yy — ^ 
at (x ,y ). Hence, we have 

u xx Uyy(x ,yo) > 0, 
a contradiction! This implies that 

SUPdU = SUPqdU. 

The same argument applies to the differential inequality: 

det(D 2 u) < 0, in D. 

(2). The Maximum principle is true for many nonlinear elliptic and 
parabolic partial differential equations of second order. It is a basic 
tool in the study of these kinds of pde's. 



72 LI MA 

9.2. MP for heat equations 1. 

Assume that D is a bounded smooth domain in R n . As before, we 

let 

Q = {(x,t);x G D,0 < t < T}, 

and let 

d p Q = {(x, t) G Q; x G dD, or t = 0} 

be the parabolic boundary of the domain Q. 
We now consider the heat equation (HE): 

Lu := u t — a 2 Au = f(x, t) 

on Q. We assume that u G C 2,1 (Q) n C(Q) is a solution to (HE). Here 
C 2,1 (Q) is the space of continuous functions which have continuous 
x-derivatives up to second order and continuous t-derivative. 

Theorem 11. Assume that u G C 2,1 (Q) fl C(Q) satisfies Lu < on 
Q. Then we have 

maxQu(x,t) = maxd p Qu(x,t). 

Proof. We first assume that Lu(x, t) < in Q. Assume that there is 
an interior point (a^o^o) G Q such that 

■u(a;o,to) = maxQu(x,t). 

Then we have at the point (xq, to), Vu = 0, D 2 u < 0, and Ut > 0. This 
implies that Au(x ,to) < and 

Lu(x ,t ) > 0, 

a contradiction. Hence, u can not have an interior maximum point and 
then 

maxQu(x,t) = maxg p Qu(x,t). 
For general case when 

Lu(x,t) < 0, (x,t) G Q, 

we let, for e > 0, 

v(x, t) = u(x, t) — et. 

Then 

Lv (x, t) = Lu(x, t) — e < 0, (x, t) G Q. 
Then we can use the result above to v(x,t) to get 

maxQv(x,t) = maxQ p Qv(x,t). 

Note that 

maxQu(x,t) < maxQv(x,t) 

and 

max dpQ v(x,t) < max dpQ u(x,t) + eT. 



EQUATIONS FROM MATHEMATICAL PHYSICS 73 

Sending e — > 0, we get 

maxQu(x,t) < maxd p Qu(x,t). 
However, since d p Q C Q, we have 

maxqu{x,t) > maxg p Qu(x,t). 

Hence, 

maxQu(x,t) = maxQ p Qu(x,t). 
We are done. □ 

Let's consider an easy application of the maximum principle. 

Theorem 12. Given f e C(<5) wi/i i 71 := maxq\f\. Assume that 
ueC 2 ' 1 (Q)f]C(Q) satisfies 

Lu = f(x, t), (x, t) e Q 

with u(x,t) = (j>(x,t) on d p Q. Then we have 

maxQu(x,t) = FT + maxQ p Q\(p(x,t)\. 

Proof. Denote by 

B = max dpQ \<f)(x,t)\. 
Consider the function w(x,t) defined by 

w(x,t) = Ft + B ±u(x,t). 

Then 

Lw = F ± / > 0, on Q 

and on d p Q, 

w>0. 

By the Maximum principle above, we have 

w(x,t) > 0, (x,t) e Q, 
which implies the conclusion wanted. □ 
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9.3. MP for heat equations 2. 

We show in this subsection that the maximum principle cab be 
proved to be true for the heat equation by another method. 

Theorem 13. Given a bounded continuous function g{x) on the smooth 
bounded domain D C R n . Let u = u(x,t) be a non-trivial smooth 
solution in C(D) to the heat equation (HE): 

u t = Au, on D x [0,T) 

with initial and boundary conditions 

u(x,0) = g(x), x G D; u(x,t) = 0, dD x [0,T) 

where D is a bounded piece-wise smooth domain in R n . Then, 

max \u(x, t)\ 

is a monotone decreasing function of t. 
Proof. Consider a new function defined by 

J k (t) = [ \ u \ 2k dx. 
Jd 

Recall that, for each t, 

M(t) := su PxeD \u(x,t)\ = lim J k (t) 1/2k . 

First, we note that 

inflimJ k (t) 1/2k < suplimJ k (t) 1/2k < M(t). 

So we only need to show that 

M(t) < inflimJ k (t) 1/2k . 

In fact, For fixed t, since a is a bounded continuous function in D, we 
have a sub-domain Q C D such that on Q , 

M(t)-S< \u{x,t)\ 

Hence 

(M(t) - 5)\Q\ 1/2k < J k (t) 1/2k < M^W™ 
Sending k — > +oo, we have 

(M(t) -S)< inflimJ k (t) 1/2k < suplimJ k (t) 1/2k < M(t). 

Since 5 > can be any small constant, we get 

lim k -+ +00 J k (t) > M(t). 

We remark that the last inequality is also true for any bounded con- 
tinuous function defined on a unbounded domain D. 
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Note that 



2k J u 2k - x u t = 2k J u^Au. 



dJk 
~dt 

By using integration by part, we have 

^ = -2k(2k - 1) J u 2k - 2 \Vu\ 2 <0. 

Hence Jk(t) is a monotone decreasing function of t. Therefore, 

max \u(x, t)\ 

xdD 

is a monotone decreasing function of t too. □ 

Remark 1. If we assume u > in Q satisfies 

u t < Au, in Q. 

The same conclusion as in the theorem above is also true. 
Remark 2. 

Given a constant A > 0. If we assume u > in Q satisfies 
Ut — a 2 Au < Au, in Q. 
Then, arguing as before, we have 

Jk(t) < Jk(0)e 2kAt , 

which implies that 

M(t) < M{0)e 2at . 

Remark 3. Note that for any bounded continuous function / defined 
on a unbounded domain D with f D \ f\ p dx < oo for some p > 0, we still 
have 

vl/fc 



M:= sup D \f(x)\ = Iim(/ \ f\ k dx)' 

k— >oo J y) 



According to the remark made in the proof of the theorem above, we 
only need to prove that 

suplim k ^{ [ \f\ k dx) 1/k < M. 
Jd 

However, it is almost trivial since 

suplimk^i [ \f\ k dx) 1/k < suplim k ^ OQ (M k - p [ \f\ p dx) 1/k = M. 
Jd Jd 

With an easy modification, we can extend the above theorem into 

unbounded domain D. We omit the detail. 
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10. LECTURE TEN 

10.1. Variational principle. 

Variational principle is a important tool in sciences. One of the 
beautiful result in the calculus of variations is the Noether principle, 
which can not be presented here. We can only touch the the concept 
of variational structure. We discuss the variational method for Laplace 
equation in a region of dimension two. 

Given a bounded domain DcR" We introduce the Dirichlet energy 
(so called the variational structure) 

I(u) = / \Vu\ 2 dx, 
Jd 

for any smooth function u : D — > R with fixed boundary value 

u(x) = <f)(x), x G 3D. 

We denote A by such a class of functions. We look for a function w E A 
such that 

I(w) = inf ueA I(u). 
Assume such a function w exists. Then for every £ G Cq(D), we have 

w + if G A 
for every t G R. So we have, for all t G R, 

/(«;) < J(iu + tg) 

and then 



-|/(«; + ^)lt=o = 0. 



Note that 



Hence, 

d 



I(w + t£) = J [\Vw\ 2 + 2tVw ■ Vf + t 2 |Vf | 2 ]. 



/(w + *0|t=o = 2 /vwVe = 0. 



Assume that w£C 2 . Then using integration by part, we have 



/ 



(A<) = 0. 



Since £ can be any testing function in Cq(D), we get the Euler-Lagrange 
equation for the functional /(•), by the variational lemma (see our 
textbook [1]), that 

Aw = 0, in D, 
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with the boundary condition w = <fi on dD. This is the famous Dirichlet 
problem in the history. The recent study is about the obstacle problem 
for harmonic functions. 

Remark: We mention that the existence of the infimum w can be 
proved by using the Hilbert space theory, which is a beautiful part of 
the course Functional Analysis. It will be really helpful for readers 
to know some knowledge about Sobolev spaces. Here are some simple 
facts about it. Let I = [a,b]. We consider L 2 (I) as the completion of 
the space C^°(I) with respect to the norm 

l/l = (^l/(^)| 2 ^) 1/2 - 

In the other word, the element / (which will be called a L 2 -function) 
in L 2 (I) is defined by a Cauchy sequence {f k } C Cq°(J), i.e., 

l/*-/il-0 

as k,j — > oo. We say that the element / G L 2 (I) has derivative 
g E L 2 (I) if there is a Cauchy sequence {fk} C Cq°(7) such that 

f k -> /, and (f k ) x -> g 

in L 2 . If this is true, we often call g the strong derivative of / and 
denote by g = f x . Then is easy to see that the integration by part 
formula 




for all 4> £ is true. People often use this formula to define the 

weak derivative of /. It can be showed that the strong derivative is the 
same weak derivative. Furthermore, one can show that both definitions 
are equivalent each other. Hence, we simply call g the derivative of /. 
We denote by Hq(I) the Sobolev space which consists of L 2 functions 
with derivatives. For /, h e Hq(I), we define 

(/, h) = J f x h x dx. 

Then one can show that (•, •) is an inner product on Hq(I), and Hq(I) 
is a Hilbert space. 
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10.2. Some comments on Nonlinear evolution problem. 

In some sense, we try to give a few comments about recent research 
directions in the study of non-linear PDE. We assume in this section 
that readers know the Banach fixed point theorem. We recall the Fixed 
Point Theorem: Let E be a Banach space and let B C E be a closed 
convex subset. Let / : B — > B be contraction mapping, i.e., there is a 
positive constant re G (0, 1) such that 

1/0*0 ~ f(y)\ < K\x-y\. 

Then there is a fixed point x* G B of /, i.e., f(x*) = x*. Here | • | is 
the norm of E. 

One often use Fixed Point Theorem to get the local existence of 
solutions to non-linear evolution problem. Let's see a little more. We 
study the evolution problem of type 

ut = Au + F(u) 

or 

u tt = Au + F(u) 

in R n with nice initial data. Using the Duhamel principle, we can 
reduce the problem into the form: 

u = S(t)(uo)+ [ S(t-r)F(u)dr := /(«). 
Jo 

Here S(t)(u ) is the solution to the problem with F — 0. Let 

5 = J B fl (5(*)M)- 
We have to show that for small R and T, the mapping / satisfies 

/ : B - B, 

i.e., is an invariant set of / and / is a contraction, that is, there 
exists a positive constant re G (0, 1) such that \f(x) — f(y)\ < re|x — y| 
for all x,y £ B. Generally speaking, to make / : B — > B, we need to 
take .E to be defined by the energy method such that E is a Banach 
algebra. Then we can use the interpolation inequality method to show 
that / is in fact a contraction mapping. 

The interesting feature in non-linear problems is that there is a blow 
up phenomenon. We show some simple examples below. 
Example One: Consider the ODE 

u t = u 2 , 

with the initial data u(0) = a > 0. It is easy to see that the solution is 
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which goes to infinity as t — > This is the blow up property of the 
ODE. 

For a non-linear evolution PDE, we can study the blow up property 
by introducing some monotone or concave function A(t) defined by the 
integration like 

pit) = minouix, t) 

or 

pit) := J pix)uix,t) k 

or 

pit) := J pix)\V x uix,t)\ k , 

where k is an integer, p(x) > is a solution of some elliptic equation. 
We can show that p(t) will satisfies some differential inequality which 
makes pit) go to negative as t — > +00. This is absurd since p(t) > 
for all t. This gives us the blow up property. 

Example Two: We consider the following nonlinear parabolic equa- 
tion 

u t = Au + u 2 , (x, t) G [0, 2tt] x [0, T) 
with the non-trivial initial data 

u(x, 0) = <f>(x) > 
and the boundary data t) = m(x + 2/r, t) for x G [0, 2tt]. Let 

pit) = min xe[ o,2n}uix,t). 

Then we have 

p\t)>pitf. 

This gives us that 

Pit) > — — 
1 — ta 

for a = min xe [o t 2ir]<f>(x) > and u has to blow up to +00 in finite time 
before aT 1 . 

Example Three: Given a constant f3 > \. We consider the following 
nonlinear parabolic equation 

u t = Au + u 2 + (3u, (x, t) G [0, 2rr] x [0, T) 

with the non-trivial initial data 

u(a;,0) = <f>{x) > 0, x G [0,2tt] 

and the boundary data w(0,t) = = u(2n,t) for x G [0, 2ir]. By the 
Maximum principle, we have 

u(x,t) > 0, x G [0,2tt] x {t > 0}. 
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Choose C — \ such that 



O I Sin(-)^. 



Let 

f 2w x 
pit) = C I sm{—)u(x,t)dx. 
Jo 2 

Then, we have 

V = C J sin(|)ut 

r-2vr „ r 2n 



C sm{^-)Au + C sin(^)« 2 + (3p(t) 

f 2n r 1 
Cjf sin(-)u 2 + ((3 - -)p(t) 



> P(t) 2 . 

This gives us the blow up as before in finite time. 

Very recently, we have used similar idea to prove the long standing 
conjecture on non-linear Schrodinger equation on compact Riemannian 
manifolds. 
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10.3. Nonlinear PDE of first order. 

This is a fun model for non-linear pde of first order. 
Let us first consider a nice example: 

u l + u l = 1, 

which is called the eikonal equation. Introduce the curve (x(t),y(t)) 
by 

dx dy 

Tt =Ux ' dt =Uy - 

Then by the equation, we have 

(^) 2 + (f) 2 = l- 

This implies that t is the arc- length parameter for the curve. Doing 
differentiation to the eikonal equation, we have 

Hence along the curve, 

d 2 x dx dy 

~d¥ = Uxx tt +Uxy dt = 

Similarly, we have 

d 2 y 

— = 0. 
dt 2 

This says that the curve is a straight line, saying, 

x(t) = at + c, y(t) =bt + d, a 2 + b 2 = 1. 
Then we have the simple equation: 

u x = a, u y = b. 

So 

u = ax + by + C, 

with a 2 + b 2 = 1. 

We now give some material for ambitious readers (who may see [2] for 
more). We consider the first order PDE of the following form (FPDE): 

^Ai(x,u)u x i = B(x,u). 

i 

We look at the graph 

z = u(x). 

Then the PDE says that the normal vector (Vu, —1) of the graph is 
orthogonal to the vector field 

(Ai(x, z), A n (x, z),B(x, z)). 
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This suggests us that the characteristic curve (x(t), z(t)) should satisfy 

If we can solve this system starting from a given surface (x(s),z(s)), 
s G r, of dimension n — 1 to get a solution x = x(s,t),z = z(s,t). 
If we solve the equation x(s,t) = x for (s = x(x),t = t(x)), then 
u = z(s(x),t(x)) is the solution to (FPDE). 

For more, please look at the great book of D.Hilbert and R. Courant 

[2]- 
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11. Appendix A 

One may consider this section as a last one to our course. We have 
three parts. One is for the variational lemma. Another is about the 
method of separation of variables. The other is about some open prob- 
lems in PDE's. Our readers may do research for these problems in the 
future. 

11.1. Variational lemma. 

Lemma 14. Let f G C°(D) for a bounded domain in R n . Assume 
that 



D 



f(x)£(x)dx = 
for every £ G Cq 3 (D). Then f = 0. 

Proof. We argue by contradiction. Assume that there is a point xq G D 
such that f(xo) ^ 0. Then we may assume that f(xo) > and there is 
a small ball B e (x ) such that f(x) ^ for every x G B e (x ). Take 

= p(\x-x \) 

where 

_ i 

p(r) = e 1 ~ r ' 2 , < r < e 
and p(r) = for all r > e. Note that p(\x - x \) G C^(D), and 
/(x)/o(|a: - x |) > 0, in B e (x ). 

Hence, 

/ /(x)e(a;)tfa = / f(x)p(\x-x \)>0. 

JD Jcg°(D) 

This gives a contradiction to our assumption of our lemma. Therefore, 
the lemma is true. □ 



Example. Assume that u = u(t) G Cq[0, L\. Define 
F{u) = j\\\ur-\u*-g{t)u{t))d^ 

where g(t) is a given continuous function on [0, L}. Then, if SF(u) = 0, 
u satisfies 

-u =u{t) 2 + g{t). 

In fact, by definition, 

O=<5F(m),0>= I {u'4>' -u 2 <P- g)dt 
Jo 
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for all (j) G C£°(0,L). Integrating by part, we get 

0= / (-«" -u 2 - g)<f)dt. 
Jo 

By variational lemma, we know that 

u" + u 2 + g = 0, 
which is equivalent to our equation above. 
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11.2. Sturm-Liouville theory. In the old time before Fourier, people 
often asked if any continuous function, which is defined on the interval 
[0, L\ with the fixed boundary condition X(0) = a,X(L) = b, can be 
expressed as the expansion of trigonometric functions. The answer is 
yes since we know the beautiful theory of Fourier series. 

In the method of separation of variables, we meet the problem if we 
can use the eigenf unctions of the following eigenvalue problem 

(p(x)X')' + [A + q(x)]X = 0, 0<x<L, 
- ai X'(0) + b l X(0) = 0), 
a 2 X'(L) + b 2 X{L) = 0, 

where the constants satisfy > 0,b{ > and + bi ^ 0, and the 
functions p(x) > and q(x) are continuous on [0, L]. 

Luckily, the answer is yes too. This is the Sturm-Liouville theory. 
Roughly speaking, the Sturm-Liouville theory says that there are an 
increasing eigenvalues \ n — * oo and eigenfunctions X n (corresponding 
to A n ) with 

L pL 

X n X m dx = 0, / \X n \ 2 dx = 1 
Jo 

for A n 7^ A„, such that every smooth function / with the boundary 
condition 

f -ai/'(0) + 6i/(0) = 0), 
1 a 2 f'(L) + b 2 f(L) = 0, 

can be expressed as 

f(x) = J2CnX n (x). 



L 
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11.3. Three open problems. 
1. De Giorgi conjecture: 

In 1978, De Giorgi proposed the following conjecture: 
If u is a solution of the scalar Ginzburg-Landau equation: 

Au + u{\ - u 2 ) = 0, on R™ 

such that \u\ < 1 and -^u > on R n , and 

lim u(x',x n ) = ±1, 

x n — »±oo 

for any x' G R n_1 , then all level set of u are hyper-planes, at least for 
n < 8. 

The case when n < 8 was solved recently by Savin, who is a PhD 
student of L. Caffarelli. For more reference, please see our paper: 

Y.H.Du and Li Ma; Some remarks related to the De Giorgi conjec- 
ture, Proc. AMS, 131(2002)2415-2422. 

2 Serrin conjecture 

In the past years, I and my student D.Z.Chen studied the following 
conjecture of Serrin: 

There is no bounded positive solution {u, v) to the following elliptic 
system on R n : 

-Am = v q , 
-Av = u p , 

where the indices p > 0, q > satisfy 

1 1 1 2 

+ > 1 - 



p+1 q + 1 n 

We have some progress in this problem. Please see our paper 
Li Ma and D.Z.Chen, A Liouville type Theorem for an integral Sys- 
tem, CPAA, 5(2006)855-859. 

3. A question from Yanyan Li and L.Nirenberg 
It is not hard to show by using the Taylor expansion that, if a non- 
negative function satisfying 

\u"(x)\ < M 

on [-R, R] , then we have 

|u'(0)| < y/2u(0)M 

for M > and 

. ,. Xl u(0) R„ 
\u (0)|<-^ + -M 



for M < 2 ^. 
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Yanyan Li and L.Nirenberg (2004) asked what is the best constant C 
in the following result in the ball Br C R n : 

Assume that u is a non-negative C 2 function in closure of the ball 
Br satisfying: 

max I Aul = M. 

Be 

Then there is a constant C depending only on n such that 

|Vu(a;)| < Cy/u(0)M 

ifi2>)/#>2|x|, 

|V«(a;)| < CC^- + RM) 
R 

\l>2\x\>R<^. 
This is an open question. 
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12. Appendix B 



This is a section for part of homework exercises in my courses in the 
past years. 

12.1. Homework 1. 

Exercise 1. Assume that the vibrating string of length L moves with 
two end points fixed at x = and x — A. Assume its initial position 
forms a triangle with the x-axis such that the middle point of the 
string is the highest point and its height is b > 0. Find the determined 
problem for this string. 

Answer: Let u be the height position of the vibrating string at 
(x,t). Then u satisfies 

u tt = a 2 u xx , for (x,t) e (0,A) x (0, oo) 

with the initial data 

u ( x > 0) = x G [0, ^] and u(x, Q) = ^(A-x), x G A] 

and the boundary data u(0,t) = and u(A,t) = 0. The relation 
between L, b, and A is 

L 2 A 2 u2 
4 4 

Exercise 2. Assume ^ x e R n , n > 2. Verify that 
A|:zf = (n/? + /?(/? - 2))|x| /3 " 2 . 

Answer: Note that d x tr = for r = \x\. 
Exercise 3. Assume < u = u(x) e C 2 . Denote by 

|V M | 2 = ^|^| 2 



(1) . Assume that / : R — > R is in C 2 . Prove that 

A/H = /'( M )A M + /"( M )|V M | 2 . 

(2) . In particular, for /3 e R, prove that 

An? = Pu^Au + (3{(3 - l)u p - 2 \Vu\ 2 . 
Exercise 4- Define the Gamma function T(s) for s > by 

POO 

r(s) = / e- r r s ' l dr. 
Jo 

Show that the definition is well-defined, and T(s + 1) = sT(s) for s > 0. 
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12.2. Homework 2. 

Exercise 1. Solve the equation 

u t + 3u x = 0, (x, t) e R x (0, oo) 

with initial data u = x 2 at t = 0. 
Answer: 

u(x, t) = (x — 3t) 2 . 
Exercise 2. Solve the equation 

ut + u x = u, (x, t) e R x (0, oo) 

with initial data u = cosx at t = 0. 
Answer: 

t) = e* cos(x — £). 
Exercise 3. Solve the equation 

2u t = — ct, (x, t) G R x (0, oo) 

with initial data u = 2xe x2 / 2 at £ = 0. 

Answer: Note that along the characteristic curve 

t 

x = ~+A, 

we have 

du . t . 

- = - xu =(-- + A)u. 

Then, 

u(x,t) = {2x + t)e x2/2 . 
Exercise 4- Solve the problem 

y Uyy x u xx 

with the initial data 

u(x, 1) = (f)(x), u y (x, 1) = ip(x). 

Answer: 

u(x,y) = -<f)(xy) + |^(-) 



+ / s 2(p(s)ds 



4 



'si/ 



2 



s 2 ip{s)ds. 



xy 
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12.3. Homework 3. 

Exercise 1. Assume that Y is a solution to 

— y" + y = 0, x > 

with the initial conditions y(0) = and y'(0) = 1. Using the Duhamel 
principle to find the expression for solution to the the following ODE 

-y" + y = g(x), x>o 

with the initial conditions y(0) = A and y'(0) = B. 
Answer: 

y(x) = AY'(x) + BY{x) - [ g(r)Y(x - r)dr. 

Jo 

Exercise 2. Given constants a±, ...,ak- Assume Y(x) is a solution to 
the the higher order ODE of the form 

y W + a 1 y« e - 1 '> + ...a k y = 

with the initial data 

1/(0) = 0, i/(0) = 0, ... 2/^-^(0) = 1. 
Find the solution expression for the higher order ODE of the form 
y {k) +a 1 y (k ' 1) + ...a k y = g(x) 

with the initial data 

y(0) = b , j/'(0) = 6 1 , ... 2 / (fc - 1) (0) = 6 fe - 1 . 
Answer: Define, for < j < k — 1, 

Uj (a;) = + aiF '" 1 )^) + ... + ajY(x). 

Then the solution expression is 

nx 

V( x ) = ^2 b j u j( x )+ g(r)Y(x - r)dr. 
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12.4. Homework 4. 

Exercise 1. Find a solution to the following Sturm-Liouville problem: 

X" + AX = 0, xe (0,L), 

and X(0) = X'(L) = 0. 
Answer: 

~ U 

and 

,(k + h)vt x 

X = 0, or X k (x) = B sin(- j^x) 

where k — 0, 1, 2 

Exercise 2. 

u tt = a 2 u xx , < x < L,t > 0, 
u{0,t) = = u x {L,t), 
u(x, 0) = sin 2 ^, Wt(x, 0) = x(L — x). 

Answer: 

an7r . .amr . nn 

u(x,t) = 2_^{A n cos(—^-t) + B n sm(—j—t))sm(—x), 

n>l 

where 

8 



n(n 2 — 4)7r 

for n — 2k + 1, and A n = for n = 2/c, and 

8L 3 



-Br; 



an% 4 



for n = 2k, and S„ = for n — 2k + 1, 
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12.5. Homework 5. 

Exercise 1. Given a domain Q C R 2 . Consider the following prob- 
lem: 

-Au = f(x,y), (x,y) G Q 

with the boundary condition u = <f> on <9f2. Find the Green function 
for the problem when 

(1) . Q is the half upper plane. 

(2) . D, = {x > 0,y > 0}. 

(3) . f2 = {— oo < x < oo,0 < y < L}, where L > is a fixed 
constant. 

Answer: 

(1) . Define 

G 1 ((x, y), (c, v )) = r((x, i/), (C, 77)) - r((x, y), (C, -77)). 

(2) . Define 

G 2 ((x, y), (C, 77)) = G l ((x, y), (C, 77)) - G 1 ((x, y), (-(, V )). 

(3) . Define 

00 

G 3 ({x f y),(C,v))= E Gi((z,J/),(C,2A;L + »7)). 
fc=— 00 

Exercise 2. Assume D = {(x,y);y > 0} is the upper half plane. 
Compute the Poisson kernel in such D. 

Answer: We just remark that on the boundary dD, the unit our 
normal is v — (0, —1). 

-k Exercise 3. Let B = B R (0) and assume G C(dB). Define 

= -— / — dcr v , 

nu n R J dB \x - y\ n 

where oo n = |-Bi(0) is the volume of the unit ball. Prove that 

Au(x) = 0, x G B. 

-k Exercise 4- Let G D be a bounded domain in R ra . Given u G 
C 2 {D) with 

Au(x) = f(x), ; x G D. 
Define D* = 1(D), where / is the inversion transform 

x ED -> /(x) = x/\x\ 2 G 

on the unit sphere {\x\ = 1}. 

(1). Prove that I 2 = identity, that is, I(I(x)) = x for x G -D. This 
implies that £> = /(£>*). 
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(2) . Let 

w(x) = \x\ 2 ~ n u(x/\x\ 2 ), x e D* 
which is called the Kelvin transform of the function u. Then 

Aw(x) = \x\ 2 ' n f(x/\x\ 2 ), x E D*. 

(3) . Using the conclusion of (2), find the Green function on B x 
Answer: (3). Define, for x,y G -Bi(O), 

G(x,y) = T(x,y)-\y\ 2 ~ n T(x,I(y)). 

Then, G(x ) y) is the Green function on -Bi(O). 
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12.6. Homework 6. 

Exercise 1. Given a constant B > 0. Compute the Fourier transform 
of the function f(x) = 1 for \x\ < B and f(x) = for \x\ > B. 



Answer: 



Exercise 2. Given two constants a > 0, c and two bounded smooth 
functions f(x,t) and Using the Fourier transform to solve the 

following problem 

u t = a 2 u xx + cu + f(x,t), — oo < x < oo, t > 0, 
0) = 

Answer: 

u(x,t) = e ct [[ K(x-U)<l>(Z)dZ 



+ l e ~ CTdT J ' K(x-£,t-T)f(Z,T)dZ]. 
We can obtain this result by two methods. 

Method One: Let w = e~ ct u and fi(x, t) = e~ ct f(x, t). Then we have 

w t = a 2 w xx + /i 

with 

w(x, t) = 4>{x). 
Method Two: Note that for u(£,t), 

u t = (-a 2 \e + c)u + f(Z,t), 

with 

-k Exercise 3. Define, for m6S, 

d 2 

H[u](x) = (-—+x 2 )u(x), 

A[u](x) = (-^- + x)u(x)i 

C[u](x) = (—^ + x)u(x). 
dx 
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Show that 

CA[u] = (H-l)[u], 

AC[u] = (H + l)[u], 

H{e~ x2 ' 2 ) = e~ x2 / 2 

H(C k (e- x ' 2/2 )) = (2k + l)C k (e- x2/2 ), k = 1,2, .... 

We remark that C k (e~ x ' are eigenfunctions of the operator H and 
the polynomials defined by 

h k (x) = e x2 ' 2 C k {e- x2 ' 2 ) 

are called the Hermite functions. H is called the Hermite operator, 
which plays an important role in modern methods of mathematical 
physics and its higher dimensional analogue 

L = -A + \x\ 2 

plays an important role in Bose-Einstein condensate, 
•k Exercise 4- Define, for u e S, a new function 

oo 

U(t) = u(t + 2im), 

n=—oo 

which is called the periodization of u. Clearly, U(t) is a periodic 
function such that U(t + 2tt) = U (t) on R. Denote by u the Fourier 
transform of u. Using Fourier series theory, you show that 



U(t) = (v^)- 1 J2 u{n)e int , 



which is called the Poisson summation formula. Then setting t — 0, 
we get 



£«(n) = (v^)- 1 £ti(n). 



•k Exercise 5. Recall that the theta function ~&(t), t > 0, is defined 

by 

oo 
n=—oo 

Using the Poisson summation, you show that 

r 1/2 0(i/f) = 

•k Exercise 6. Assume that 
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Find B(x). 

Answer: First we write 

e 



-s(i+\e) ds 



Then we have 



Recall that 



i + i«i 2 

B{x)= j^S?'' ds L eii ^ 



e iay-sy* dy =^l/2 e -£ 
R s 



Then we have 



e i*<s\e d t. = (![)n/2 e -M 



L 

Hence, we have 

B{x)=2- n ' 2 / s^e"'"^^, 

which is called the second Bessel potential. 
* Exercise 7. Given (3 > 0. Assume that 

MO = (1 + |^|2^/2 • 

Find Bp{x). 

Hint: One can find the expression of Bp in E.M. Stein's book, 
Princeton Univ. press,1970, or W. P.Ziemer's book, GTM120, Springer, 1989. 
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12.7. Homework 7. 

Exercise 1. Let D be a bounded domain in R n . Consider u G 
C 2 (D) n C(D) such that 

-Au = f(x), x G D 
u{x) = <f)(x), x G 3D. 

Then 

max^|w(x)| < maxaz)|0(x)| + — sup|/(a;)|, 
where d = diam(D) := sup x ^ & qd\x — y\. 
Proof. We denote 

F — sup\f(x)\, M = max dD \<f)(x)\ 

Take p G D. Let 
Note that 

£(x)±u(x) > 

on 3D and 

-A(£(x)± «(:*;)) >0 
Applying the Maximum principle to the function 

w(x) = £(x) ±u(x), 

we get that w(x) > on D. This implies that 



maxjj\u(x)\ < maxdD\(j)(x)\ + — sup\f(x)\. 

Zi it 



□ 
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12.8. Homework 8. 

Exercise 1. Let Q = {(x, t); < x < L, < t < T}. Let a > be a 
constant. Assume that u G C 2 ' 1 ^) satisfies 

u t = a 2 u xx + f(x, t), < x < L,0 < t < T, 
u (0,t) =0 =u x (L,t), 
u(x, 0) = 

Show that maxq\u t {x, t)\ can be bounded by the suitable norm of / 
and 0. 

Answer: There is a uniform constant C — C{T, a) such that 
maxQ\u t \ < C[\f\ c i { Q } + \<j)"\ C [o,L]}- 

Exercise 2. Given <fi G C(D) and a(x, t) G C(Q) with A = maxq\a(x, t)\. 
Assume that u G C 2,1 (Q) n C(<5) satisfies 

Lu := u t — a 2 Au = —u 2 + a(x, t)u, (x, t) G Q 

and the initial data 

u(x,0) = (j)(x) > 0, 
and the boundary data u(x, t) = for x G 3D. Prove that on Q, 

< u(x,t) < Msup D \<p(x)\ 

for some constant M > depending only on n, T, and A. 

Proof. Assume that u < somewhere, saying, 

u(xi,ti) = —(eti + e 2 ) 2 

for some arbitrary small e and ti > 0. Consider 

w = u + et + e 2 . 

At (xi, ti), we have 

ti) = 0. 

Note that 

Aw = Au 

and at t = 0, w = <j) + e 2 > 0. Hence, there is a small < t < ti 
such that w(x,t) > for (x,t)D x (0,to] and for some minimum point 
x G -D for u> such that 

w(x ,to) = 0, Au(x ,t ) = Aw(x ,t ) > 0, 

which implies that 

u(x ,t ) = -(et + e 2 ). 

Hence, we have 

Lw(x ,t ) < 0, 
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which gives us that 

Lu(x ,t ) < -e, 
Then, using the equation for u at (xo,to), 

-e > (et + e 2 ) 2 - a(x ,t )(et + e 2 ). 

However, this is impossible for small e and t\. Hence, we have 

u > 0, on Q. 

Once we have w > 0, we can use the maximum principle 2 (see 
Section 9.3) to conclude that 

u(x,t) < M supd4>{x). 

□ 

Remark. One may try another method below. Assume that m : = 
u(x 2 ,t 2 ) = maxqU > 2sup£)(p. Then at (x 2 ,t 2 ), we have 

Lu > 0. 

Using the equation, we get 

— m 2 + a(x 2 ,t 2 )m > 0. 

Hence, we have 

m 2 < Am 

and m < A. So this method can not be used to give the result wanted. 
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12.9. Homework 9. 

Exercise 1. Let F and G be two smooth functions in R. Assume 
that u(x,t) E C^R x (0, T)) is a bounded function such that 

lim x ^ ±00 u(x,t) = 

uniformly in t and 

d t F(u) + d x G(u) = 0. 
Here d x G(u) = G'(u)u x . Show that 

— / F(w)d:r = 0. 
dtJ R K J 

We call the quantity f R F(u)dx the conservation law. 

Exercise 2. Let G be a vector-valued smooth function in R. Assume 
that u(x,t) G C^R™ x (0, T)) is a bounded function such that 

Zim^i^oo^a;, i) = 

uniformly in t and 

d t u + div x G(u) = 0. 

Show that 

d 
dt 

for each k = 1, 2, 3.... 

^ Exercise 3. Using the conclusions of Exercise 1 and Exercise 2 to 
find at least 5 conservation laws of the KdV equation: 

u t + uu x + u xxx = 0. 

Hint: Define 

E n (u)= [ {[u^] 2 + c n u[u^] 2 -q n {u,...M x n - 2) )dx. 
Jr 

for suitable constants c n and polynomials q n . Here denotes the 
x-derivative of order n for the function u. 

* Exercise 4- Let u = u(x,t) be a smooth solution to the equation 

u t - u xx + cu 2 x = 0. 

Find a function such that the function w = (f>(u) satisfies the heat 
equation 

w t - w xx = 0. 

Answer: w = e~ cu , which is called the Cole-Hopf transformation, 
•k Exercise 5. Consider the Burgers' equation 

Ut ~ u xx + uu x = 0, x e R 



/ u k dx = 

JR" 
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with initial data u = f(x). Let 



Show that w satisfies 



/X 
u(y,t)dy. 
-co 



w t - w xx + ^w 2 x = 



and use the conclusion of Exercise 4 to find an explicit expression of 
u. 

Answer: 

r , s -\x-y\ 2 /o(«) , 

v(r f) _ k( x ~y) e 4t _ 2 dy 

a \ X l L ) - -l x - vl 2 fQ(y) 

tj R e « 2 dy 

where / (x) = f* f(y)dy. 
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12.10. Homework 10. 

Exercise 1. As an exercise for pleasure, we invite readers to find the 
Euler- Lagrange equation for the area functional 



A(u) 



[ v 7 ! + \Vu\ 2 dx 
Jd 



among the class A above. 

Answer: From this variational structure, we know the famous 
minimal surface equation: 

-div( , VU =) = 0. 

There is a rich literature about this equation. 

Exercise 2. Try to use the method of separation variables to the 
Minimal surface equation in dimension two in the following way. As- 
sume u = u(x, y) = f(x) + g(y) satisfies 

d x ( . ^ =) + 8y{ . ^ =) = 0. 

Then we have 

/ 9 



0. 



i + f( x y i + g'( y y 

Prove that for some constant A > 0, 

1 cos(Ay) 

uixu) = —log\ t—^-\. 

v 1 A y| cos(Aa;) 1 

This solution is called the Sherk's surface. 

Exercise 3. Find solution u(x,y) of the form u = f(r), r 
\J x 2 + y 2 , to the equation 



u 



Answer: 

f(y) = jCosh{Ar + B). 
Exercise 4- Find the Euler-Lagrange equation for the integral 



Answer: 



I(u) — / L(x, u, u x i, u x n)dx. 
Jd 



d u L - ^2 l d xi {d Uxj L) = 0. 
j'=i 
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* Exercise 5. Find the Euler-Lagrange equation for the integral 

Vu\ 2 dx 



Jtu) = ! 
Jd 

on the space Cq(D) with the constraint 

/ |-u| 2 <ir = 1 
Jd 

Answer: 



—Am = Am. 
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Summary 

This is a brief lecture note for the courses on " Introduction to Equa- 
tions from Mathematical Physics" or "Methods for Mathematical Physics" . 
Our motivation to write this note is to show the mainstream methods 
used in solving some typical partial differential equations like wave 
equation, heat equations, and the Laplace equation. 

For students who would like to take the final examination of this 
course, please look my homepage at 
http : / /faculty . math . tsinghua . edu . cn/~lma 

for the highlight of the course. 

Thank you very much for your patience in attending 
my courses. 
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